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Abstract 

We introduce a nonperturbative interaction for light-cone fluctuations containing 
quarks and gluons. The qq interaction squeezes the transverse size of these fluctuations 
in the photon and one does not need to simulate this effect via effective quark masses. 
The strength of this interaction is fixed by data. Data on diffractive dissociation of 
hadrons and photons show that the nonperturbative interaction of gluons is much 
stronger. We fix the parameters for the nonperturbative quark-gluon interaction by 
data for diffractive dissociation to large masses (triple- Pomeron regime). This allows 
us to predict nuclear shadowing for gluons which turns out to be not as strong as 
perturbative QCD predicts. We expect a delayed onset of gluon shadowing at x < 10~^ 
shadowing of quarks. Gluon shadowing turns out to be nearly scale invariant up 
to virtualities ~ AGeV^ due to presence of a semihard scale characterizing the 
strong nonperturbative interaction of gluons. We use the same concept to improve 
our description of gluon bremsstrahlung which is related to the distribution function 
for a quark-gluon fluctuation and the interaction cross section of a qqG fluctuation with 
a nucleon. We expect the nonperturbative interaction to suppress dramatically the 
gluon radiation at small transverse momenta compared to perturbative calculations. 



1 Introduction 



The hght-cone representation introduced in [|I|] is nowadays a popular and powerful tool to 
study the dynamics of photo-induced (real and virtual) reactions. The central concept of 
this approach is the non-normalized distribution amplitude of qq fluctuations of the photon 
in the mixed (p, a) representation, where p is the transverse qq separation and a is the 
fraction of the light-cone momentum of the photon carried by the quark (antiquark). For 
transversely and longitudinally polarized photons it reads 0] , 



^l\p,c.) = ^xO^^'xMep) (1) 

Here x ^"^^ X are the spinors of the quark and antiquark respectively. Ko{ep) is the modified 
Bessel function, where 

= a(l - a)Q^ + m\ . (2) 

This is a generalization of 0] to the case of virtual photons [|], 
The operators O^'^ have the form, 

= (J ■ e + z(l - 2a) (a ■ n) (e ■ Vp) + (a x e) ■ Vp , (3) 

6^ = 2Qa{\- a)a -n , (4) 

where the dimension- two operator Vp acts on the transverse coordinate p; n = p/pis a. unit 
vector parallel to the photon momentum; e is the polarization vector of the photon. 

The advantage of the light-cone approach is the factorized form of the interaction cross 
section which is given by the sum of the cross sections for different fluctuations weighted 
by the probabilities of these Fock states [1], |, |]. The flavor independent color-dipole cross 
section aqq first introduced in ||^ as dependent only on transverse qq separation p. It vanishes 
quadratically at p ^ due to color screening. 



= C(p,3)p^ (5) 



where C(p, s) is a smooth function of separation and energy. In fact, C(p, s) also depends on 
relative sharing by the q and q of the total light cone momentum. We drop this dependence 
in what follows unless it is important {e.g. for diffractive gluon radiation). It was first 
evaluated assuming no energy dependence in pQCD 0, ^ and phenomenologically ^ at 
medium large energies and p's and turned out to be C ~ 3. There are several models for 
the function C(p, s) [e.g. in |TU|, |rT|), unfortunately neither seems to be reliable. In this 
paper we concentrate on the principal problems how to include nonperturbative effects, and 
do not try to optimize the form of the cross section. For practical applications it can be 
corrected as soon as a more reliable model for C (p, s) is available. We modify one of the 
models mentioned above |TD[ which keeps the calculations simple to make it more realistic 



and use it throughout this paper. 
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The distribution amplitudes (|I|) control the mean transverse qq separation in a virtual 
photon, 

(p^) ^- (6) 

Thus, even a highly virtual photon can create a large size qq fluctuation with large prob- 
ability provided that a (or 1 — a) is very small, aQ"^ ~ m^. This observation is central 
to the aligned jet model p[. At small soft hadronic fluctuations become dominant at 
any a. In this case the perturbative distribution functions which are based on several 
assumptions including asymptotic freedom, are irrelevant. One should expect that nonper- 
turbative interactions modify (squeeze) the distribution of transverse separations of the qq 
pair. In Section 2.1 we introduce a nonperturbative interaction between the quark and an- 
tiquark into the Schrodinger type equation for the Green function of the qq pair |]12|, |T3|, 
The shape of the real part of this potential is adjusted to reproduce the light-cone wave 
function of the p-meson. We derive new light-cone distribution functions for the interacting 
qq fluctuations of a photon, which coincide with the known perturbative ones in the limit 
of vanishing interaction. The strength of the nonperturbative interaction can be fixed by 
comparison with data sensitive to the transverse size of the fluctuations. The observables 
we have chosen in Section 2.2 are the total photoabsorption cross sections on protons and 
nuclei and the cross section for diffractive dissociation of a photon into a qq pair. 

For gluon bremsstrahlung we expect the transverse separation in a quark-gluon fluctua- 
tion to be of the order of the typical color correlation length ~ 0.3 fm obtained by several 
QCD analyses |]T5[ - |T^. This corresponds to the radius of a constituent quark in many 
effective models. To the extend that the typical q — G separation is smaller than the qq one 
we expect gluon radiation to be suppressed. This results in particular in a suppression of 
diffractive gluon radiation, i.e. of the triple- Pomeron coupling, which is seen indeed in the 
data. 

In Section 3.1 we assume a similar shape for the quark-gluon potential as for the qq one, 
but with different parameters. A new light-cone distribution function for a quark-gluon 
fluctuation of a quark is derived, which correctly reproduces the known limit of perturbative 
QCD. 

Comparison with data on diffractive excitations with large mass fixes the strength of the 
nonperturbative interaction of gluons. An intuitive physical picture of diffraction, as well 
as a simple calculation of the cross sections of different diffractive reactions is presented in 
Appendix A. A more formal treatment of the same diffractive reactions via calculation of 
Feynman diagrams is described in Appendix B. 

A crude estimate of the interaction parameters is given in Section 3.1.1 within the 
additive quark model (AQM). For this purpose the cross section of diffractive gluon radiation 
by a quark, qN — > qGN, is calculated in Appendices A. 2 and B.l, based on general 
properties of diffraction (Appendix A.l) and the direct calculation of Feynman diagrams. 

Quite a substantial deviation from the results for the AQM is found in Section 3.1.2 and 
Appendix C where the diffractive excitation of a nucleon via gluon radiation, N N ^ X N 
is calculated. The high precision of the data for this reaction allows to fix the strength of 
the nonperturbative interaction of gluons rather precisely. 
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The cross sections of diffractive gluon radiation by mesons and photons are calculated 
in Appendices A. 3 and B.2. In Section 3.1.3 we compare the values of the triple-Pomeron 
couplings (calculated in Appendix C) for diffractive dissociation of a photon and different 
hadrons and find a violation of Regge factorization by about a factor of two. 

Our results for the cross section of diffractive dissociation 7* qqG N in the limit of 

vanishing nonperturbative interaction can be compared with previous perturbative calcula- 
tions [|18|, |19|. In this limit we are in agreement with |]19|, but disagree with [^.Q The source 
of error in [^] is the application of Eq. (|A.6|) to an exclusive channel and a renormalization 
recipe based on a probabilistic treatment of diffraction. 

Diffractive radiation of photons is considered in Appendices A. 4 and B.3. It is shown 
that no radiation occurs without transverse momentum transfer to the quark (in contrast 
to gluon radiation). Therefore, the cross section for diffractive production of Drell-Yan 
pairs is suppressed compared to the expectation of |^ which is also based on an improper 
application of Eq. ( |A.6| ) to an exclusive channel. 

Section 3.2 is devoted to nuclear shadowing for the gluon distribution function at small x. 
Calculations for many hard reactions on nuclei (DIS, high pt jets, heavy flavor production, 
etc.) desperately need the gluon distribution function for nuclei which is expected to be 
shadowed at small x. Many approaches |21|-|31| to predict nuclear shadowing for gluons can 



be found in the literature (see recent review ^^). Our approach is based on Gribov's theory 
of inelastic shadowing |3^ and is close to that in 13^, |31| which utilizes the results [|^, ^ 
for the gluonic component of the diffractive structure function assuming factorization and 
using available data. Instead, we fix the parameters of the nonperturbative interaction using 
data on diffraction of protons and real photons. Besides, we achieved substantial progress in 
understanding the evolution of diffractively produced intermediate states in nuclear matter. 

Nuclear suppression of the gluon density which looks like a result of gluon fusion G G — > 
G in the infinite momentum frame of the nucleus, should be interpreted as usual nuclear 
shadowing for the total interaction cross section of fluctuations containing a gluon if seen in 
the rest frame of the nucleus. We perform calculations for longitudinally polarized photons 
which are known to be a good probe for the gluon distribution function. Although the 
physics of nuclear shadowing and diffraction are closely related, even a good knowledge 
of single diffractive cross section and mass distributions is not sufficient to predict nuclear 
shadowing completely, but only the lowest order shadowing correction. A technique for 
inclusion of the multiple scattering corrections was developed in |T2|, [13] which includes 



evolution of the intermediate states propagating through the nucleus. These corrections 
are especially important for gluon shadowing which does not saturate even at very small 
X in contrast to shadowing of quarks. In Section 3.2.1 we find quite a steep x-dependence 
of gluon shadowing at > AGeV'^ which is rather weak compared to what have been 
estimated in |31[]. Shadowing starts at smaller values of x < 0.01 compared to the 
shadowing of quarks. Such a delayed onset of gluon shadowing is a result of enlarged mass 
of the fluctuations containing gluons. 

As soon as our approach incorporates the nonperturbative effects we are in position 



^In spite of the claim in that their result coincides with that of [ij], they are quite different. We are 
thankful to Mark Wiisthoff for discussion of this controversy. 



4 



to calculate shadowing for soft gluons as well. This is done in Section 3.2.2 using two 
methods. In hadronic basis one can relate the shadowing term in the total hadron-nucleus 
cross section to the known diffractive dissociation cross section. This also give the scale for 
the effective absorption cross section. A better way is to apply the Green function approach 
which includes the nonperturbative gluon interaction fixed by comparison with data for 
diffraction. With both methods we have arrived at a similar shadowing, but the Green 
function approach leads to a delayed onset of shadowing starting at x < 0.01. We conclude 
that gluon shadowing is nearly scale independent up to ~ AGeV"^. 

The nonperturbative interaction of the radiated gluons especially affects their transverse 
momentum distribution. One can expect a substantial suppression of radiation with small 

related to large transverse separations in quark-gluon fluctuations of the projectile quark. 
Indeed, in Section 3.3 we have found suppression by almost two orders of magnitude for 
radiation at /ct = compared to the perturbative QCD predictions. The difference remains 
quite large up to a few GeV of momentum transfer. Especially strong nonperturbative 
effects we expect for the /cr-distribution of gluon bremsstrahlung by a quark propagating 
through a nucleus. Instead of a sharp peak at A;t = predicted by pQCD now we 
expect a minimum. 

2 Virtual photoproduction of quark pairs 

2.1 Green function of an interacting quark-ant iquark pair 

Propagation within a medium of an interacting qq pair which has been produced with 
initial separation p = from a virtual photon at a point with longitudinal coordinate Zi 
and developed a separation p at the point Z2 (see Fig. [l|) can be described by a light-cone 
Green function Gqq{zi,pi = 0; Z2, P2 = p)- The evolution equation for this Green function 




Figure 1: Illustration for the Green function Gqq{zi, pi = 0;Z2,P2 = p) 
for an interacting qq fluctuation of a photon, as defined by Eq. (0). 
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was studied in 0-1111], 



i-^Gqq{zi,pi]Z2,p2) 



A. 



2pa{l 



a] 



+ Vqq{z2j,a) 



Gqq{zi, pi] Z2, P2) ■ 



(7) 



The first term on the r.h.s. is analogous to the kinetic term in a Schrodinger equation. It 
takes care of the phase shift for the propagating qq pair. Indeed, the relevant phase factor 
is given by exp[i dz qL{z)], with the relative longitudinal momentum transfer q^. The 
latter is defined by 



''^'^ = 2-p = 2pail - a) 

Here p is the photon momentum; M is the effective mass of the qq pair (which varies with 
z) and is the photon virtuality. It depends on the transverse momentum fc-r of the quark 
(antiquark) which is replaced by the Laplacian, — Ap, in the coordinate representation 

The imaginary part of the potential Vqq{z2, p,(y) is responsible for absorption in the 
medium which is supposed to be cold nuclear matter. 



lmVqq{z2,p,a) 



(T. 



qq 



(P) 



Pa{Z2} 



(9) 



Here Pa{z) is the nuclear density and we omit the dependence on the nuclear impact pa- 
rameter. (Tqq{p, s) is the total interaction cross section of a colorless qq pair with a nucleon 
introduce in (|^). Eq. (|^) with the imaginary potential was used in |]12[ to calculate 
nuclear shadowing in deep-inelastic scattering. In other applications the quarks were treated 
as free, what is justified only in the domain of validity of perturbative QCD. 

Our objective here is to include explicitly the nonperturbative interaction between the 
quarks in (|^. We are going to rely on a nonrelativistic potential, which, however, should 
be modified to be a function of the light-cone variables p and a. This general problem is, 
however, not yet solved. Nevertheless, we try to model the real part of the potential based 
on its general properties. Particularly, the qq pair is supposed to have bound states which 
are vector mesons. 

It is assumed usually that the wave function of a vector meson in the ground state 
depends on p and a according to 



^'y(p,a) = f{a) exp 



(10) 



In order for this to be a solution of (0) the real part of the potential should be, 

ReVgq{z2,p,a) 



2pa{l — a) 



^our Green function is related to that in jlj] by 
Gqq{zi,pi = 0;z2,p2 = p) = exp[~ie'^{z2 - zi)/2pa{l - a)] W{zi,pi = 0; Z2,p2 = p) 
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Unfortunately, no reliable way to fix the form of a{a) is known. A parameterization pop 
ular in the literature is a(a) 



2 a \Ja{l — a), which results from attempts to construct a 
relativistic approach to the problem of a qq bound state ( see and references therein). 
In this case, however, the mean qq separation p (x \ j \J ail — a) increases unrestrictedly 
towards the endpoints a = 0, 1. Such a behavior contradicts the concept of confinement 
and should be corrected. The simplest way to do so is to add a constant term to a{a) 
(the real form of a{a) may be quite different, but so far data allow only for a simple two 
parameter fit). 



a [a] = Oq + 4a^ a(l 



One can roughly evaluate by demanding that even at a 
does not exceed the confinement radius, 



a) . (12) 
0, 1 the transverse qq separation 



-1 



Aqcd , (13) 



i.e. ao ~ 200 Me Comparison with data (see below) leads to a somewhat smaller value. 

In what follows we study the consequences of the interaction between q and q in the form 
dO) ^ (HH) quark wave function of the photon, and we discuss several observables. 

Let us denote the Green function of a qq pair propagating in vacuum (ImV^ = 0) as 



Gqq{zi, Pi] ^2, f>2)- The solution of (|^) has the form ||3^ 



Gqq{Zl, Pi] Z2, P2) 



2 7C i sinfu; Az] 
i e 



exp 



i a^(a) 
sinfcu Az) 



(Pi + pI) cos(u; Az) - 2 pi ■ p2 



^2 Az 



2 p a{l — a) 



where Az = Z2 — zi and 



a) 



(14) 



(15) 



p a{l 

The normalization factor here is fixed by the condition Gqq{zi, pi] Z2, P2)\z2=zi = S^{pi — f>2)- 
Now we are in the position to calculate the distribution function of a qq fluctuation of a 
photon including the interaction. It is given by the integral of the Green function over the 
longitudinal coordinate zi of the point at which the photon forms the qq pair (see Fig. 1), 



C/(p,«) 



iZq^e 



QQ 



4:7ipa{l — a) 



dZl (X Gqq 



Zl, Pi] ^2, P2/ 



Pi=0; P2=p 



(16) 



The operators O^'^ are defined in (§)-(^. Here they act on the coordinate pi. 
If we write the transverse part as 

xd^X = A + B-Vp, , 

then the distribution functions read, 

^^^(p,a) = Zqy/^ [A<l>o(e,p,A) + i?$i(e,p,A) 



(17) 



;i8) 
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where 



A = A(a) 



The functions $0,1 in (P^-(P^) are defined as 



'^'o(e,P,A) 



47r 



dt 



sh(At) 



exp 



2^2 



Stt 



00 


A 


2 


/* 




exp 


_sh(At)_ 



4 

AeV 
4 



cth(At) - 1 



cth(At) - 1 



(19) 
(20) 

(21) 
(22) 



Note that the q — q interaction emerges in ([T8|)-(p^ through the parameter A defined 
in (pOl). In the hmit A — {i.e. — >■ 0, a is fixed, a 7^ or 1) we get the well known 
perturbative expressions (|I]) for the distribution functions, 



^o(e,P,A) 
<^i(e,P,A)) 



A=o 27r 
ep 



Ki{ep) 



1 



A=o 2-np ^ ' 2tx 
In contrast to these relations, in the general case, i.e. for A 7^ 

$i(e,p,A) ^-V$o(e,p,A) . 



Viro(ep) 



(23) 
(24) 

(25) 



In the strong interaction limit A ^ iriq (or if both (Q^, mq — > 0)) which is appropriate 
particularly for real photons and massless quarks, the functions $0,1 acquire again simple 
analytical forms, 

$o(e,P,A) 



1 



A^oo 47r 



-a [a)p 



$i(e,p,A) 



P 



■ exp 



The interaction confines even massless quarks within a finite range of p. 



(26) 
(27) 



2.2 Absorption cross section for virtual photons 

For highly virtual photons, ^ a^(a), according to (|2Dp A ^ and the effects related to 
the nonperturbative qq interaction should be gone. Although for very asymmetric configu- 
rations, a{l — a) -C 1, see (||) the transverse qq separation increases and one may expect 
the nonperturbative interaction to be at work, it does not happen if the dipole cross section 
is independent of p at large p. 

Thus, our equations show a smooth transition between the formalism of perturbative 
QCD valid at high and our model for low where nonperturbative effects are important. 
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The absorption cross sections for transversely (T) and longitudinally (L) polarized virtual 
photons , including the nonperturbative effects read, 



a 



tot 



i 

2N,Y,Zlaemj da J d^pagg{p,s) jm^ $2(e, p, A) + [a' + (1 - a)^] |$i 



.p,\) 



tot 



1 

SQ"^ NcJ2^q'^em J da a {1 - a) J d'^pagg{p, s) 



P,A) 



(28) 
(29) 



Here Nc is the number of colors, and the contributions of different flavors F are summed 
up. 

According to (||) the dipole cross section vanishes <Jqq{p, s) oc at small p <^ 1 fm. Such 
a behavior approximately describes e.g. the observed hierarchy of hadronic cross sections as 
functions of the mean hadronic radii pj. We expect, however, that the dipole cross section 
flattens off at larger separations p > 1 fm. Therefore, the approximation crgg{p) oc p^ is 
quite crude for the large separations typical for soft reactions. Even the simple two-gluon 
approximation [^, ^ provides only a logarithmic growth at large p f^, and confinement 
implies a cross section which becomes constant at large p. Besides, the energy dependence 
of the dipole cross section is stronger at small p than at large p [^. We use hereafter a 



parameterization similar to one suggested in 10 



1 — exp 



(30) 



where po{s) = 0.88 /m (sq/s)"'^^ and Sq = 1000 GeV^^. In contrast to [0] all values depend 
on energy (as it is supposed to be) rather than on x and we introduce an energy dependent 
parameter cro(s), 

5plis)\ 



tot\ 



1 + 



(31) 



42 



otherwise one fails to reproduce hadronic cross sections. Here {r'^h)-^ = 0.44 ± 0.01 fm^ 
is the mean square of the pion charge radius. Cross section (^) averaged with the pion 
wave function squared automatically reproduces the pion-proton cross section. We use the 
results of the fit for the Pomeron part of the cross section. 



(^tot{s) 



23.6 (s/so)''-"' mb 



(32) 



where sq = 1000 Ge^^. We fixed the parameters comparing data with the proton structure 
function calculated using Eqs. (p8D-(p9|) and the cross section (^). Agreement is reasonably 
good up to ~ 20 GeV"^ sufficient for our purposes. 

To fix from data the parameters ao,i of the potential we concentrate on real photoab- 
sorption which is most sensitive to nonperturbative corrections. The photoabsorption cross 
section with free quark fluctuations in the photon diverges logarithmically at rUg — > 0, 



T 

^tot - 



CToln 



nig po, 



(33) 
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Inclusion of interaction between the quarks in the photon makes the photoabsorption 



cross section finite at 



0. 



tot 



0(Xi) - (\)[X2) 



(34) 



where 



and 



1 + QqPo 

2 2 ' 



X2 



2 ' 



X 



(x) = 4 In ( - ) - 2 x + (4 + x) vT+x In 



'^l + x + 1 



(35) 



(36) 



In this case the cross section of photoabsorption is independent of the quark mass in the 
limit rriq/ao,! <^ 1. 

We adjust the values of Oq and ai to the value of the photoabsorption cross section 
aZ = IGOpb at ^/s = 200 GeV |g, ||]. Eq. (|3|) alone does not allow to fix the two 
parameters ao and ai completely, but it provides a relation between them. We found a 
simple way to parameterize this ambiguity. If we choose 



,1.15 



112 MeV f 



:i-v- 



1.15 



lQ5MeVf 



(37) 



the total photoabsorption cross section, turns out to be constant (within 1%) if v varies 
between and 1. This covers all possible choices for oq and ai. 

In order to fix v in (^) one needs additional experimental information. We have tried 
a comparison with the following data: 



(i) The cross section of forward diffraction dissociation —>■ qqN (the PPR graph in 
the triple- Reggeon phenomenology ||46[|), 



da{'jN qqN) 



dt 



t=o 



1671 



da d^p 



a\p) . 



(ii) The total photoabsorption cross sections for nuclei (high-energy limit). 



(^tot 



2 J d'^b J dad^p '^gg{a,p) 



exp 



l^cr{p)T{b) 



where 



T{b) 



dzpAib,z) 



(38) 



(39) 



(40) 



is the nuclear thickness function and the nuclear density pA{b, z) depends on impact param- 
eter b and longitudinal coordinate z. This expression can be used for virtual photons as well 
with a proper discrimination between transverse and longitudinal photons. 
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A calculation of the observables (i) and (ii) shows, however, a surprising stability of the 
results against variation of v in (pT]): the cross sections change only within ~ 1% if f varies 
between and 1. Thus, we were unable to constrain the parameters Oq and Oi any further. 

We have also calculated the effective interaction cross section of a qq pair with a nucleon, 



qqN 
eff 



J da (fp 


^qq{a,p) 


2 

a2(p, s) 


_ (^^) 


J da d^p 


^qq{a,p) 


2 

(t(p, s) 





(41) 



which is usually used to characterize shadowing for the interaction of the qq fluctuation of a 
real photon with a nucleus {e.g. see in ^]). We got at o"f//^ = 30 m6 at ^/s = 200 GeV. 
This well corresponds to the pion- nucleon cross section (|3^) 0"^? = 31.7 mb at this energy. 
This result might be treated as success of VDM. On the other hand, A calculation [E7| using 



VDM and 0"^* ~ 25 mb instead of afoi at lower energy for photoproduction of p-mesons off 
nuclei is in good agreement with recent HERMES measurements [^]. 

However, a word of caution is in order. The nucleus to nucleon ratio of total photoab- 
sorption cross sections in the approximation of frozen fluctuations (reasonably good at very 
small x) reads p, [3^, [31[| , 



a 



7*A 
tot 



An 



7*Af 
tot 



2 

V) 



d^b ( 1 — exp 



m 



(42) 



Expanding the exponential up to the next term after the double scattering one (1/4) cie// 
one gets (1/24) (cr^)/(cr). This is 1.5 times larger than (1/24) (Tgjj/(cr) if to use the dipole 
approximation a (x. (? and a Gaussian distribution over p for color triplet (g — q) or color 
octet [G — qq) dipoles. 



3 Gluon bremsstrahlung 

3.1 Radiation of interacting gluons 

In processes with radiation of gluons, like 

q^N g + G + X 
7* + A^ ^ g + g + G + X 



(43) 
(44) 



the interaction between the radiated gluon and the parent quark traveling in nearly the same 
direction may be important and significantly change the radiation cross section and the 
transverse momentum distribution compared to perturbative QCD calculations 0, pl[] . 

We describe the differential cross section of gluon radiation in a quark-nucleon collision 
in the factorized light-cone approach [BH] 



(iV(g -^qG) _ 1 
rf(lna) d^kT ~ M 



rfVirfV2exp[i/ET(ri-f2)]^^g(a,fi)^'Gg(a,r2)aG(ri,f2,a) , (45) 
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where 

crG(n,^2,a) = ^{(^Gqq{ri,ri-ar2)+aGqq{r2,r2-ari)-agg[a{fi-f2)]-aGG^^^ ■ (46) 

Hereafter we assume all cross sections to depend on energy, but do not show it explicitly 
for the sake of brevity (unless it is important). 

The cross section of a colorless Gqq system with a nucleon (JGqq{ri,r2) is expressed in 
terms of the usual qq dipole cross sections, 

9 r ~| 1 

(^Gqqin, ra) = -{aqqin) + cTqg{r2)\-- agq{fi - fa) (47) 

fi and 7^2 are the transverse separations gluon - quark and gluon - antiquark respectively. 
In ( ^61) (TGGii^) = I c"5g('^) is the total cross section of a colorless GG dipole with a nucleon. 
The cross sections of reactions (|4^ ) - (|4^) integrated over have simple form. 



da(q — >• qG) 
d{\na) 




da{'y* qqG) 



d{\n.a 



G) 



1 




= ^ doLg J d^R 




aG<l 





X 



d'r 



^5G(i? + f, ac) Vgg(^ + + 



qG[r,aG) 



'gg( 



- Re v^;g(^ «g) ^qGiR + r, a^) <(i? + r) + <(r) - a'^^iR) 



(49) 



Here is the fraction of the quark momentum carried by the gluon; R and r are the 
quark-antiquark and quark-gluon transverse separations respectively. The three terms in 
the curly brackets in (^) correspond to the radiation of the gluon by the quark, by the 
antiquark and to their interference respectively. 

The key ingredient of (|^), (|48|) and (149|) is the distribution function \1/Gg(tt,r) of the 
quark-gluon fluctuation, where a is the fraction of the light-cone momentum of the parent 
quark carried by the gluon, and f is the transverse quark-gluon separation. This function 
has a form |53[ similar to (|I|), 



free 



(50) 



where the operator T is defined in [pq 



[fix a) + a e* ■ (a X V) - i{2 - a) e* ■ V 



(51) 



We treat the gluons as massless since we incorporate the nonperturbative interaction ex- 
plicitly and do not need to introduce any effective mass. 
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The factor r differs from e as defined in (Q), 

r2 = a'^ml (52) 

In the general case the distribution function including the interaction between the quark 
and gluon can be found via the Green function Ggdzi, pi, Z2, P2) for the propagation of a 
quark-gluon pair, in analogy to (p!6D, 



i J as/ 3 y /- -"T \ 

^gG(p,a) = 7^ -jz V / dziixT x)GqGizi,pi;z2,p2) ^ _ ^ (53) 

2TTpa{l — aj J ^ pi=0; P2=p 

—00 

Let us add a few comments as to why this direct analogy holds. Eq. (^) might give the 
impression that we would have to implement the interaction between all three partons: the 
gluon, the quark and the antiquark. Checking the way in which this equation was derived, 
one realizes, however, that this is not the case. We studied gluon bremsstrahlung from a 
single quark and then expressed the radiation amplitude as a difference between the inelastic 
amplitudes for a qG system and an individual q. This is how acqq has to be interpreted and 
this is why one should only take the q — G nonperturbative interaction into account. 

The evolution equation for the Green function of an interacting qG pair originating 
from the parent quark at the point with longitudinal coordinate zi with initial transverse 
separation pi = looks similar to (|^ with the replacement e ^ r and V^q(z2,P, «) ^ 
V^g(^2,P, «)• We parameterize the quark-gluon potential in the same way as in (0) for 
quark-antiquark. 

Re V,g{z2, p, a) = / ^ , , (54) 
2pa[l — a) 

where 6^(a) = fog + 4&i a (1 - a). 

The solution of the evolution equation for the quark-gluon Green function in absence of 
absorption (IrnVqc = 0) looks the same as ( p!^ with replacement a{a) ^ b{a). 

The following transformations go along with (0) - (pTj). The vertex function in (^) is 
represented as, 

X Fx = ^ + ^ • Vp, , (55) 



then the result of integration in (|53|) is 



vl/,G(p,a)=2^(^^j [D$o(r,p,A) + E-$i(r,p,A)] , (56) 

The functions $o(t, p, A) and $i(r, p, A) are defined in ( pTf ) - (|2^). However, A is now defined 
by 

A^^. (57) 

One might argue that the quark-gluon potential we need (and which we shortly shall 
constrain by comparison with experimental data) could simply be obtained by adding two 
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quark-quark potentials with an appropriate color factor. Such a procedure could, however, 
lead to a completely wrong results as we want to illustrate by the following example. 

Motivated by perturbative QCD one might expect that the gluon-gluon and quark-quark 
potentials differ simply by a factor 9/4 (the ratio of the Casimir factors). However, this 
relation is affected by non-trivial properties of the QCD vacuum which makes the interaction 
of gluons much stronger [U, The octet string tension Kg is related to Op, the slope of 
the Pomeron trajectory in the same way as the color triplet string tension relates to the 
slope of the meson Regge trajectories |l55| . 



^^f, = ,^-^GeV/fm. (58) 

ZTTttp 

Here a'p = 0.25 GeV~^. Thus, the value of Ks is in fact four times larger than the well 
known qq string tension, ^3 = 1 GeV/ fm, and not only by a factor 9/4. 

Another piece of information about the strength of the gluon interaction which supports 
this observation comes from data on diffractive dissociation. The triple-Pomeron coupling 
turns out to be rather small |56]. If interpreted as a product of the Pomeron flux times 



the Pomeron-proton total cross section, the latter turns out to be an order of magnitude 
smaller than the proton-proton one. Naively one would again assume that the Pomeron as 
a colorless gluonic dipole should interact 9/4 times stronger that an analogous qq dipole 
(such a consideration led some authors to the conclusion that gluons are shadowed at small 
X in nuclei stronger than sea quarks). The only way to explain this discrepancy is to assume 
that the gluon-gluon dipole is much smaller. This, in turn, demands a stronger gluon-gluon 
interaction. Thus, diffraction is sensitive to nonperturbative interaction of gluons. We shall 
use this observation in in the next section to fix the corresponding parameters 60 and 61. 

3.1.1 Diffractive bremsstrahlung of gluons. The triple-Pomeron coupling in 
the additive quark model 

Let us start with diffractive dissociation of a quark, qN — > qGN. We assume the diffractive 
amplitude to be proportional to the gluon density G{x,Q'^) = xgi^XjQ"^) [|56|, |5^ (it should 
be a non-diagonal distribution if the energy is not very large) as is shown in Fig. ^ Since the 
amplitude is predominantly imaginary at high energies one can use the generalized unitarity 



relation known as Cutkosky rule [58 



2lmAab = Y.^acAl , (59) 

c 

where Aab is the amplitude of the process a ^ b and a, b denote all the particles in initial 
and final states respectively. In the case under discussion a = {q, N}, b = {q, G, N}, and c 
denotes either ci = {g, A^|} or C2 = {q,G,N^}, where is a color-octet excitation of the 
nucleon resulting from gluon radiation(absorption) by a nucleon. 

In what follows, we concentrate on forward diffraction, i.e. the transverse momentum 
transfer qr = 0. In this case the diffractive amplitude reads, 

72 

F(a, kT,qT = 0) = ^l^ e^'^^^^^cia, p) a{p) , (60) 
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where /ct is the transverse momentum of the radiated gluon and 

9 

^(P) = g^5g(p) ■ (61) 

Eq. ( |60D is derived in Appendix A. 2 in a simple and intuitive way based on the general 
properties of a diffractive process discussed in Appendix A.l. A more formal derivation 
based on a direct calculation of Feynman diagrams and the Cutkosky rule ( ^9] ) is presented 
in Appendix B.l. 

The relation ( [BDD is valid for any value of a. In contrast to the inclusive cross section for 
gluon bremsstrahlung the diffractive cross section depends on a only via the distribution 
function. 
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The amplitude (|60|) is normalized according to 

da{qN qGN) 



d{lna) (Pkx (Pqx 



F{a,kT,qT) 



The distribution for the effective mass squared, 
has the form, 

M2 



T 



a) 



da{qN qGN) 



dM"^ dql 

which transforms in the limit 

M2 da{qN qGN) 



27r^ 



'?T=0 

oo into 



dJxirp 



F{a,kT,qT) 



(62) 

k'^/a, at qr = 



(63) 



dM'^ dql 



(JT=0 



lim 

167r a^O 



d^p 



'^qG{a,p)(r{p, s) 



(64) 



where s = Mg/^, where Mq = 1 Gel^^ and ^ = xp = 1 — xp ^ M'^/s. 

Since dissociation into large mass states is dominated by the triple-Pomeron (3P) graph 
the value on the l.h.s. of (0) is the effective 3P coupling G^p^qN — > XN) (see definition in 
46| ) at qp = 0. It can be evaluated using G^p{pp — > Xp) ^ l.bmb/GeV'^ as it follows from 



measurement by the CDF Collaboration |59| (according to |6T| we divided the value of G^p 
given in |^ by factor 2). This value is twice as small as one derived from the triple- Regge 
analyses at medium large energies. This is supposed to be due to absorptive corrections 
which grossly diminish the survival probability of large rapidity gaps at high energies [3C]. 
One can see energy dependence of Gap even in the energy range of the CDF experiment 

m. 



Assuming the additive quark model (AQM) to be valid one can write. 



r^AQM 
^3P 



{qN XN) ^ - GspiNN XN) ^ 0.5 
3 



mb 



(65) 



GeV^ 

(see below about interference effects) . To compare with this estimate we calculate the triple- 
Pomeron coupling 
section (p 



4|) using the distribution function in the form (|56| ) and the dipole cross 
Gf^^iqN 



XN) 



27 a. 



8nJ 



In 



(66) 



where t\ = 6^(0)/2, ^2 = ^1 + 1/Po ^3 ~ 2^1^2/(1^1 + h)- The parameters (Jq and po are 
defined in (|30|). We use here a fixed value of = 0.6 which is an appropriate approximation 
for a soft process. 

Comparison of Eq. (|66|) with the value (|65D leads to a rough evaluation of the parameter 
60 of our potential (we are not sensitive to bi since keep a small). 



6^Q*^(0) ^ 570 Mel/ 



(67) 



Thus, a typical quark-gluon separation is ~ 1/^(0) ~ 0.4 /m what is roughly the radius 
of a 'constituent' quark. Note that a substantial modification of ( |65D by interference of 
radiation amplitudes for different quarks is possible. 

■^We thank Doug Jansen and Thomas Nunnemann who helped to clarify this point. 
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3.1.2 Diffractive excitation of nucleons, N N — > X N, beyond the AQM. 

The amplitude of diffractive gluon radiation N N —>■ 3qG N can be represented as a super- 
position of radiation by different quarks as shown in Fig. ^ In this process the colorless 




Figure 3: Contributions from projectile valence quarks to the amplitude 
of diffractive gluon emission in NN collisions. Six additional grapfis 
resulting from the permutation {1^2} and {1^3} have not been 
plotted. 



3q system (|3g)i) converts into a color-octet final state (|3g)8). There are two independent 
octet |3g) states which differ from each other by their symmetry under a permutation of the 
color indices of the quarks. Correspondingly, the amplitude for the process A^A^ — > \3q)sGN 
is a superposition of two amplitudes (see below). 

The contribution to the amplitude of the first graph in Fig. reads, 



F\NN 3qGN) 



abc 



V3 



((3g)8|ri'^ l(3g)i) Hpi) $({r, a}, pi, aa) 



(68) 



The second and third graphs in Fig. ^ give correspondingly. 



F'\NN 3qGN) 



abc 



V3 



((3g)8|rf rf |(3g)i) 



and 



F''\NN 3qGN) 



abc 



0"(P2) - cr(n - ^2) 



(^{ps) - - ^3) 



$({r,a},pi,aG') 



(69) 



(70) 

Here {r, a} = (ri, r2, "^3; aia2, as); are the positions of the quarks in the impact parameter 
plane; pi = p — r^, where p is the position of the gluon; 



^{{f,a},pi,aG/ai) = ^iv_,3g({r, a}) ^Gg(Pi, "g) ; 



(71) 



fabc is the structure constant of the color group, where "a" is the color index of the radiated 
gluon, and we sum over "6" and "c" . The Gell-Mann matrices A* = 2 r* act on the color 
index of i-th quark. 
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Using the relation, 

fabc rj;'^ ri^) |3g)i = {rj^^ ri^) + rf r^^) \?>q), 
one can present the sum of the amphtudes F^, F^^ and F^^^ in the form, 
F^^\NN 3qGN) = F^ + F^^ + F^^^ 



(72) 



abc 



V3 



(73) 



where Sj^ = a{p — fi) + a{p — rj)—a{fi — fj). The index "1" in F*^-'^^(A^A^ 3qGN) indicates 
that the gluon is radiated by the quark qi in accordance with Fig. |[ 

The amphtudes F^'^^ and F^^^ is obviously related to F^^^ by replacement 1^2,3. Note 
that the color structure fate 

rm ^(3) 

which is not present in ([T3|) is not independent due to 

the relation, 

fa,c (r(^) r(2) + r(2) r^^ + r^'^ r^^)) |3g)i = . (74) 

Thus, we are left with only two independent color structures, as was mentioned above. 

The full amplitude for diffractive gluon radiation squared \F{NN ?)qGN)\'^ = \F^^^ + 
^(2) _|_ ir(3)|2^ summed over all color states of the 3qG system, reads. 



E 



F{NN 3qGN) 



.i=l 



where 



A«({f},p) 
i?(i'^)({f},p) 



^12 + ^13 + ^12^13 , 



2 ^12 + Sl2(S 



13 T ^23 j 



13^23 



(75) 



(76) 
(77) 



The expressions are obtained by simply changing the indices. 

The effective triple-Pomeron coupling results from integrating ( [75D over phase space. 



GspiNN NX) = — [ d^ri d^r gc/Vs d'^p dai da2 da^ 6(fi + f 2 + fa) 

IGvr J 



X 5(1 — ai — a2 — as) ^ 

/ 



F{NN ^ 3qGN) 



(7J 



To evaluate G3p{NN NX) we use ( |19|) for "^gcip, a) and a Gaussian parameterization 
for the valence quark distribution in the nucleon. 



l"^ N^3qi{r, a})f oc exp 



u=l 



(79) 
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where {r'^h)p ~ 0.79±0.03 fm^ is the mean square radius of the proton [^|. At this point one 
has to introduce some specific model for the distributions. Quite some proposals can be 
found in the literature, and a quantitative analysis will require careful numerical studies. For 
a first qualitative discussion we make the simple ansatz for the quark momentum distribution 
in the nucleon, F^{ai, 02,0^3) oc n^(f^i ~" 1/3) which allows to continue our calculations 

i 

analytically. The details of the integration of (^) can be found in Appendix C. G^p is 
a function of the parameter 6(0). As a trial value we choose Eq. (|67D , 6(0) = 570 MeV, 
(estimated using the result of the additive quark model G^p^qN -h> XN) ^ 0.5 mb/GeV'^) 
we arrive at G3p{NN — > XN) ^ 2Amb/GeV'^ . This value is substantially higher than the 
experimental value G3p{NN XN) = 1.5mh/GeV'^. This is an obvious manifestation 
of simplifying approximations (the quark additivity) we have done. In order to fit the 
experimental value of G3p{NN —>■ NX) after the contributions of the second and third 
graphs in Fig. |^ are included we should use in (|78D 



6(0) = 650 MeV . (80) 

With this value ( |66[ ) gives, 

mh 1 

G^piqN XN) ^ 0.3 ^ - G^piNN ^ XN) , (81) 

which shows a substantial deviation from the AQM. 



3.1.3 Diffractive gluon radiation by a (virtual) photon and mesons. Breakdown 
of Regge factorization 

One can use a similar technique to calculate the cross section for diffractive gluon radiation 
by a photon and mesons. The diffraction amplitude 7(M) — qqG N is described by the 
four diagrams depicted in Fig. ^. The first two diagrams correspond to the AQM. In this 
approximation the forward (gr = 0) amplitude 7 A^ — ^ qqG N reads. 



F{q N ^ qGN) - F(qN ^ qG N) 



qG 



a 



(82) 



where = p — ri (i=l,2). p, ri 2 are the radius-vectors of the gluon, quark and antiquark 
respectively. The limit — > is assumed. 

After addition of the last two graphs in Fig. ^ the amplitude takes the form (we do not 
write out its trivial color structure). 



F{-fN ^ qqGN) 



X 



^qq{a,pi-p2) ^gcf— ,Pl 

a{pi) + a{p2) - o{pi - p2 



^ ^ — a 



(83) 
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Figure 4: Diagrams for the diffractive radiation of a gluon in photon- 
nucleon interaction, '~f*N qqGN . 



The detailed calculation of the diagrams depicted in Fig. ^ is presented in Appendix B.2. 
A much simpler and more intuitive derivation of Eq. (|83| ) is suggested in Appendix A. 3. 

If one neglects the nonperturbative effects in Eq. ( p3D {b{a) = 0) this expression coincides 
with Eq. (3.4) in |T^, but is quite different from the cross section of diffractive gluon 
radiation derived in |TB[ (Eq. (60)) (see footnote^). A crucial step in |TB| is the transition 
from Fock states which are the eigenstates of interaction, to the physical state basis. Such 
a rotation of the S'-matrix leads to a renormalization of the probability amplitudes for the 



Fock states (see Appendix A.l), rather than just the probabilities as it was assumed in 
The amplitude (|83|) is normalized as, 



da 



dilnac) dqj 



9t=0 



/ d^pid^P2da 



F{-fN qqGN) 



Direct comparison of of the cross section for diffractive gluon radiation by a photon calcu- 
lated with this expression with data is complicated by contribution of diffraction to qq states 
and nondiffractive (Reggeon) mechanisms. This is why one should first perform a detailed 
triple-Regge analysis of data and then to compare (|8^) with the effective triple-Pomeron 
coupling. Good data for photon diffraction are available at lab. energy = 100 GeV [Q. 
At this energy, however, there is no true triple-Regge region which demands s/M"^ ^ 1 
and ^ IGeV^. Therefore the results of the triple-Regge analysis in ^3| cannot be 



trusted. It is much more appropriate to use available data from HERA, particularly those 
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in IQ at y/s = 200 GeV where a triple- Regge analysis taking into account four graphs was 
performed. The result for the effective triple-Pomeron coupling 

GYp{0) = (8.19 ± 1.6 ± 1.34 ± 2.22) j^b/GeV^ (85) 

should be compared with our prediction Glp{0) = 9 fib/GeV^. To estimate the mean 
energy for the dipole cross section s/M"^ GeV^ we used the mid value = 250 GeV^ of 
the interval of measured in which corresponds to xp = 0.0064. Thus, high-energy 
data for gluon radiation in diffractive dissociation of protons and photons give the value 
(|80|) for the parameter of nonperturbative quark-gluon interaction. 

Note that the relative role of "additive" (# 1,2 in Fig. ^ and # 1 in Fig. |^) and "non- 
additive" (t^ 3,4 in Fig. ^ and 7^ 2,3 in Fig. ^ graphs depends on the relation between the 
three characteristic sizes Rh = ^ {rfj), Po and 1/6(0). In the limit Rh ^ Po, 1/^(0) the 
contribution of the "non-additive" graphs vanishes and the additive quark model becomes 
a good approximation. However, at realistic values of ~ 1 fm the "additive" and "non- 
additive" contributions are of the same order and the latter becomes dominant for small 
Rh- Particularly, this explains why the factorization relation, 

AMh N^XN) = ^''^^^^Zn^^ = Const , (86) 

atot{hN) 

i.e. independent of h, is substantially broken. We expect, 

A3p{NN^XN) = 0.025 Gel/-^ , 

A3p{ttN^XN) = 0.031 GeV-\ 

AspiKN^XN) = 0.042 Ge1/-2, 

A3p{-fN^XN) = 0.052 GeV-\ (87) 

We see that our predictions for the triple-Pomeron vertex as defined from diffractive disso- 
ciation of nucleons and photons are different by almost factor of three. On top of that, the 
absorptive corrections which are known to be larger for diffraction than for elastic scattering 
also contribute to the breaking of Regge factorization. A manifestation of these correction 
shows up as deviation between the data and the Regge based expectations for the energy 
dependence of the diffractive cross section |65 |. 



3.2 Gluon shadowing in nuclei 



It is known since long time that the parton distribution in nuclei is shadowed at small x 
due to parton fusion. In QCD this effect corresponds to the nonlinear term in the evolution 
equation responsible for gluon recombination [^T], This phenomenon is very important 
as soon as one calculates the cross section of a hard reaction (gluon radiation with high kr, 
prompt photons, Drell-Yan reaction, heavy flavor production, etc.) assuming factorization. 
Nuclear shadowing of sea quarks is well measured in DIS, but for gluons it is poorly known. 
One desperately needs to know it to provide predictions for the high-energy nuclear colliders, 
RHIC and LHC. 



21 



The interpretation of nuclear shadowing depends on the choice of the reference frame. 
In the infinite momentum frame of the nucleus it looks like parton fusion. Indeed, the 
longitudinal spread of the valence quarks in the bound nucleons, as well as the internucleon 
distances, are subject to Lorentz contraction. Therefore the nucleons are spatially well 
separated. However, the longitudinal spread of partons at small x contracts much less 
because they have an x times smaller Lorentz factor. Therefore, such partons can overlap 



and fuse even if they originate from different nucleons ||6^. Fusion of two gluons into a 
qq pair leads to shadowing of sea quarks. If two gluons fuse to a single gluon it results in 
shadowing of gluons. 

The same phenomenon looks quite differently in the rest frame of the nucleus, as shad- 
owing of long-living hadronic fluctuations of the virtual photon. This resembles the ordinary 
nuclear shadowing for the total cross sections of hadron-nucleus interaction. Indeed, the to- 
tal virtual photoabsorption cross section is proportional to the structure function F2{x, Q^). 
However, one can calculate in this way only shadowing of quarks. To predict shadowing 



of gluons it was suggested in |Q to replace the photon by a hypothetical particle probing 
gluons. Assuming for the GG fluctuation of this particle the same distribution function 
as for qq one may conclude that the effective absorption cross section providing shadowing 
is 9/4 times larger than for a qq fluctuation of a photon. Such a simple result cannot be 
true because of the strong gluon-gluon interaction which makes their distribution function 
quite different ("squeezes" it). Besides, the spin structure of the GG distribution function 
is different too. 

3.2.1 Nuclear shadowing for longitudinal photons 

Longitudinally polarized photons are known to be a good probe for the gluon structure 
function. Indeed, the aligned jet model cannot be applied in this case since the distri- 
bution function for longitudinal photons (|l|), (§) suppresses the asymmetric qq fluctuations 
with a —>■ 0, 1. Therefore, the transverse separation of the qq pair is small ~ l/Q^ and 
nuclear shadowing can be only due to shadowing of gluons. One can also see that from the 
expression for the cross section of a small size dipole ^Tj , 



4'^(rT,x) ^^as{Q')GM{x,Q') , (88) 

where Gn{x,Q'^) = xg{x,Q^) is the gluon density and ~ I/'^t- Thus, we expect nearly 
the same nuclear shadowing at large for the longitudinal photoabsorption cross section 
and for the gluon distribution, 

^kix,Q') ^ Ga{x,Q') 

aUx,Q') ^ Gn{x,Q^) ^ ' 

The estimate for nuclear shadowing for longitudinally polarized photons follows. 

Nuclear shadowing for photons corresponds to the inelastic nuclear shadowing as it 
was introduced for hadrons by Gribov 30 years ago Therefore, the term Aa{'~f*A) = 



atot{l*A) — Aa{'y*N) representing shadowing in the total photoabsorption cross section is 
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proportional to the diffractive dissociation cross section 7* ^ X [p3| , |69|1 , considered 
above. In the lowest order in the intensity of XN interaction the shadowing correction 
reads, 



d^a{YN XN) 



dMj^ dq^ 



9t=0 



X dzi dz2Q{z2 - zi) pAib,Zi) pAib,Z2)exp -iqL{z2 - zi) 



(90) 



where 



2v 



(91) 



Here v is the photon energy; Z\ and Z2 are the longitudinal coordinates of the nucleons 
Xx and X21 respectively, participating in the diffractive transition 7* X^ ^ X Xx and back 

XX2^1*X2. 

The longitudinal momentum transfer (|9ll) controls the lifetime (coherence time of 
the hadronic fluctuation of the photon, = ^/ql- It is known only if the mass matrix is 
diagonal, i.e. the fluctuations have definite masses. However, in this case the interaction 
cross section of the fluctuation has no definite value. Then one faces a problem of calculation 
of nuclear attenuation for the intermediate state X via interaction with the nuclear medium. 

This problem can be settled using the Green function formalism developed above in 
Section 2.1 |jl2|, |14[. One should switch to the quark-gluon representation for the produced 
state X = \qq), \qqG), , \qq2G), .... As one can see below an exact solution is not an easy 
problem even for the two lowest Fock states. For higher states containing two or more gluons 
it may be solved in the double-leading-log approximation which neglects the size of the 
previous Fock state and treats a multi-gluon fluctuation as a color octet-octet dipole. This 
is actually what we do in what follows, except the Fock state with only one gluon leads to a 
1/M^ mass distribution for diffraction, while inclusion of multi-gluon components makes it 
slightly steeper. This is not a big effect, besides, the nuclear formfactor substantially cuts off 
the reachable mass interval (see below). Therefore, we restrict the following consideration 
by the first two Fock states. 

For the lowest state \qq) one can write. 



SvrRe / dM]^ 



d^a{YN XX) 



dMi dq^ 



exp 



- Re y d^kx J da 




{kr, a) 



qT=0 



exp 



-iqL {Z2 - Zx, 



^ ^ k\ 

{Z2 - Zx) 



2 u ail — a) 



^Rey dVi d^r2 J da F^*^gq(r2, a) G^g(f2, ^2; n, zi) 



7*— MjgV'"!; ^) ; 



(92) 



where e was defined in (||). 
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The amphtudes of diffraction 7* ^ X in the transverse momentum and coordinate 
representations are related by Fourier transform, 



{kr, a) 



1 

2^ 



This amplitude in the coordinate representation has a factorized form, 



'^qqir,a)agg{r) . 



(93) 



(94) 



G^g(r2, Z2]ri, Zi) in (^) is the Green function of a free propagation of the qq pair between 



points Zi and Z2. It is a solution of Eq. (|^) without interaction. 



G^qi^ h r L, Zi] 



(2vr) 



(i^kr exp 



-i A;r ■ ("^2 - n) 



1 A;|, {z2 - zi) 

2 z/afl — a) 



The boundary condition for the Green function is. 



Gi(r2,2;2;ri,zi; 



qq\ 



5{r2 - ri) 



(95) 



(96) 



In Eq. (95 the phase shift on the distance Z2 — zi is controlled by the transverse momen- 
tum squared as one could expect from Eqs. (P(]|)-(|^) where it depends on the M'j^- However, 
Eq. (j9^ is written now in the coordinate representation and contains no uncertainty with 
the absorption cross section, as different from (ppj). In order to include the effects of ab- 
sorption of the intermediate state X into (p2D one should replace the free Green function 



^55(^2, 2:2; ri, 2ri) by the solution of the Schrodinger equation (^ with imaginary potential 
(|^). This was done in paper and the results have demonstrated a substantial deviation 
of nuclear shadowing from usually used approximations for transverse photons. One should 
also include the real part of the potential which takes into account the nonperturbative in- 
teraction between q and q as it is described in Section 2.1. This is important only for nuclear 
shadowing of transverse photons and and low longitudinal photons and is beyond the 
scopes of present paper. Therefore, we skip further discussion of nuclear shadowing for the 
\qq) pair and switch to the next Fock component \qqG). 

For the intermediate state ( p| ) X = qqG Eq. (0) is modified as. 



d^a{YN ^ XN) 



J ^ dMl dq^ 

J d^y2 d'^xi d^yi J dag dln^ac) 



cos[qL {Z2 - Zi)] 



qT=0 



X Fl 



.qqGi^2, m, ttg, "g) GggG^Xi, ^2, Z2] Xi, jji, Zi) F^*-.ggG(^l, Vl, "g, "g) , (97) 



where Ug and ac are the fractions of the photon light cone momentum carried by the quark 
and gluon respectively. The amplitude of diffraction 7* X X N depends on the q-q 
transverse separation x and the distance y from the gluon to the center of gravity of the qq 
pair (we switch to these variables from the previously used pi 2 for the sake of convenience, 
it simplifies the expression for kinetic energy). 
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The Schrodinger equation for the Green function Gggc describing propagation of the qqG 
system through a medium including interaction with the environment as well as between 
the constituent has the form, 

i GqqGix2, m, Z2; xi, yi, Zl) 

= 1^ - ^ A(x2) - ^ r A(f2) + V{x2,y2,Z2,aq,aG)\ 

X ^5,0(^2,^2, 2:2; ^1,^1,^1) , (98) 
with the boundary condition, 

GqqG{.X2,y2,Z2]Xi,yi,zi) _ = 6{x2 - xi) 5{y2 - yi) . (99) 

Z2 — 2^1 



The imaginary part of the potential V{x2,y2, Z2, aq, ac) in (0) is proportional to the 
interaction cross section for the qqG system with a nucleon. 



2lmV{x2,y2,Z2,aq,aG) 
1 9 



(-^q -A ( -. (^q 



y - x \ +aqq[y- X 



PA{b,z). (100) 



The real part of this potential responsible for the nonperturbative interaction between the 
quarks and gluon is discussed below. 

If the potential lmV{x2, y2, Z2, aq, ac) is a bi-linear function of x and y then Eq. (^ can 
be solved analytically. Nevertheless, the general case of nuclear shadowing for a three-parton 
system is quite complicated and we should simplify the problem. 

Let us consider nuclear shadowing for longitudinally polarized photons with high Q^. 
The latter means that one can neglect the eikonal attenuation for the qq Fock component 
of the longitudinal photon, i.e. 

:^AG{Ta) ^lGeV\ (101) 

where C is the factor in Eq. (^ and (Ta) is the mean nuclear thickness function. 

As different from the case of transversely polarized photons which distribution function 
(0) " (0) contains qq pairs with large separation (a 0, 1) even at large Q^, in longitudinally 
polarized photons small size (~ l/Q) qq pairs always dominate |l|, ||. This property suggest 
a few simplifications for the following calculations. 

1. One can neglect at large the nonperturbative qq interaction and use the perturba- 
tive photon wave function (|^) - (^). 

2. One can simplify the expression for the diffractive amplitude 7* —>■ qqGN intro- 
duced in ( P7| ) relying on smallness of the typical qq separation \x\ ~ l/Q in comparison 
with the distance between the qq and the gluon \y ~ I/60 ~ 0.3 /m. 
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3. One can also simplify the equation ( pS]) for the Green function Gqqc fixing x = in the 
expression ( p9D for the nonperturbative potential ImK(x2, ^2, ^2, ctg, «g)- This leads 



in (|98|) to a factorized dependence on variables x and y. 
As a result of these approximations and ac ^ we arrive at, 



where 



and 



(^Gcyr, s) 



9 



(102) 

(103) 
(104) 



As soon as we neglect the size of the color-octet qq pair, it interacts a gluon, this is why 
one can replace aqqG by the dipole cross section o"gg- The latter is larger than aqq by the 
Casimir factor 9/4. 

In this case the tree-body Green function factorizes to a product of two-body ones. 



GqqG{X2, m, Z2] Xl, ^1, Zl) Gqq{X2, Z2] Xl, Zi) GGG^m, Z2] Vl, Zl] 



(105) 



where Gqq{x2, Z2; xi, zi) is the "free" Green function of the qq pair, and GGG{y2, Z2',yi, zi) 
describes propagation of the GG dipole which constituents interact with each other, as well 
as with the nuclear medium. 
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dZ2 



GGG{y2,z2;yi,zi] 



2uaG(l - ac, 



+ Vim, Z2) 



GGG{y2,Z2;yi,zi) , 



where 



2lmV{y,z) = -acGiy) PA^b, z) 



On analogy to (|TT]) we assume the real part of the potential has a form 

h^y^ 



ReV{y,z) 



2 z/ 00(1 - "g) 



(106) 
(107) 

(108) 



where b Bq. 

To simplify the estimate we assume that o"gg('^! ■s) ~ Ggg{s) where Cgg{s) 
dcTGGir, s)/dr^ 



r=0- 



The solution of Eq. ( |106D has a form, 

GGGiy2,Z2;yi,zi) 
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(109) 
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where 



A 
Q 
Az 



¥ - i aoil - ao) Cqg Pa 
iA 



aG(l - olg) V 

Z2- Zi . 



The quark-gluon wave function in ( 102 ) has a form similar to 

62 
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(110) 



(111) 



Now we have all the components of the amplitude ( |102| ) which we need to calculate the 
nuclear shadowing correction (p7D. Integration in xi_2 and yi^2 can be performed analytically. 
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(112) 
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u 
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i X rriN Az ; 
A 



tcosh{VLAz) +sinh(riA^) ; 

(1 + t^) sinh(^] Az) + 2t cosh(fi A^) 



(113) 



The rest integration in (|1 1 2| ) can be performed numerically. We calculated the ratio 
-^A/AT — GA{x,Q'^)/aGN{x,Q^) for the gluon distribution functions for small values of 
Bjorken lO"'' < x < 10~^ and high = lOGeV^^ We found R^/n almost independent of 

at higher Q^. The results are depicted in Fig. ||. 

One can see that in contrast to the quark distribution the onset of nuclear shadowing for 
gluons starts at quite small x ~ lO"^. This is because the photon fluctuations containing 
gluons are heavier than qq fluctuations. Correspondingly, the lifetime of such fluctuations 
is shorter (or qi is smaller) and they need a smaller x to expose coherent effects like nuclear 
shadowing. One can expect an antishadowing effect at medium x ~ 0.1 like in F2{x,Q'^) 
which should push the crossing point Ga{x,Q'^)/G]^{x,Q'^) = 1 down to smaller x. Discus- 
sion of the dynamics of antishadowing (swelling of bound nucleons, etc.) goes beyond the 
scopes of this paper. 



27 




X 



Figure 5: Ratio of the gluon distribution functions in nuclei (carbon, 
copper and lead) and nucleons at small Bjorken x and = AGeV'^ 
(solid curves) and iOGeV"^ (dashed curves). 



A similar approach to the problem of gluon shadowing is developed in which re- 
lates shadowing to the diffractive radiation of gluons. Note that a delayed onset of gluon 
shadowing (at x < 0.02) is also expected in [^. However, this is a result of an ad hoc param- 



eterization for antishadowing, rather than calculations. The phase shift factor cos{qLAz) 
which controls the onset of shadowing in (pT]), ( |112|) is neglected in |^ assuming that 
X is sufficiently small. However, nuclear shadowing for gluons does not saturate even at 
very small x because of the form of the mass dependence of diffractive radiation of 

gluons (triple-Pomeron diffraction). The smaller the x = Q'^ /2mMV is, the higher masses 
are allowed by the nuclear form factor {q^ = (Q^ + M^)/2z/ < 1/Ra) to contribute to the 
shadowing. 

Our results also show that R'^/j^ steeply decreases down to small x and seems to have 
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a tendency to become negative. It would not be surprising for heavy nuclei if our shadow- 
ing correction corresponded to double scattering term only. However, the expression ( p^ ) 
includes all the higher order rescattering terms. The source of the trouble is the obvious 
breaking down of the unitarity limit Cdifj < Utot- This problem is well known and easily 
fixed by introducing the unitarity or absorptive corrections which substantially slow down 
the growth of the diffractive cross section. Available data for diffraction pp — >■ pX clearly 
demonstrate the effect of unitarity corrections One may expect that at very high 



energies the relative fraction of diffraction decreases. We restrict ourselves with this word 
of caution in the present paper and postpone a further study of unitarity effects for a sepa- 
rate publication, as well as the effects of higher Fock components containing more than one 
gluon. Those corrections also become more important at small x. 

Note that quite a strong nuclear suppression for gluons G^{x,Q^)/G'^{x,Q'^) < 
F^{x,Q^)/F^ {XyQ"^) was predicted in basing on the fact that the cross section of a 
color octet-octet dipole contains the factor 9/4 compared to cr^g. However, it is argued above 
in section 3.1 and confirmed by the following calculations that the observed smallness of 
the diffractive cross section of gluon radiation shows that that the strong nonperturbative 
interaction of gluons substantially reduces the size of fluctuations including the gluon. The 
situation is much more complicated and cannot be reduced to a simple factor 9/4. 

A perspective method for calculation of nuclear shadowing for gluons was suggested in 
the recent publication . Experimental data for diffractive charm production can be used 
to estimate the effect. This seems to be more reliable than pure theoretical calculations 
performed above. Indeed, the transverse separation of a heavy flavored QQ pair is small 
even at low Q^, and may be assumed to be much smaller that the mean distance between 
the QQ and the gluon. Unfortunately, the available data obtained at HERA have quite 
poor accuracy. The results from HI |71[ and ZEUS experiments are different by almost 



factor of two. Besides, the theoretical analysis @, ^ which is needed to reconstruct the 
diffractive cross section of charm production from production of D* in a limited phase space, 
introduces substantial uncertainty. According to the realistic solutions for the diffractive 



charm production differ by a factor of flve. In this circumstances we suppose our calculation 
for nuclear shadowing of gluons seems to be more reliable. 

Note that we expect much weaker nuclear shadowing for gluons than it was predicted in 
23, m g. For instance at x = lO'^ and Q^ = AGeV^ we expect Ga/AGn ^ 0.9, while a 
much stronger suppression Ga/AGn ~ 0.6 even Ga/AGn ~ 0.3 ^ was predicted 
for A ^ 200 at Q^ = AGeV^. 

It is instructive to compare the gluon shadowing at high Q^ with what one expects 
for hadronic reactions at much smaller virtualities. One should expect more shadowing at 
smaller Q^, however, the soft gluon shadowing evaluated in the next section turns out to be 
much weaker than one predicted in ^ at high Q^. 

At the same time, quite a different approach to the problem of gluon shadowing based 
on the nonlinear GLR evolution equation used in led to the results pretty close to 



ours. 
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3.2.2 Nuclear shadowing for soft gluons 

(i) Hadronic diffraction and gluon sfiadowing 

The hadron-nucleus total cross section is known to be subject to usual Glauber (eikonal) [[73l 
shadowing and Gribov's inelastic corrections |3^. Those corrections are controlled by the 



cross section of diffractive dissociation of the projectile hadron h N ^ X N which contains 
particularly the triple-Pomeron contribution. The latter as was shown above is related to 
gluon shadowing in nuclei. Namely, absorption of the incoming hadron can be treated as a 
result of interaction with the gluon cloud (in the infinite momentum frame of the nucleus) 
of bound nucleons at small x. A substantial part of this absorption is reproduced by the 
eikonal approximation which assumes the gluon density to be proportional to the number 
of bound nucleons. However, evolution of the gluon density including gluon fusion (see [EB 



and I^T], |6^ for high Q ) results in reduction of the gluon density compared to one used 



in the eikonal approximation. Such a reduction makes nuclear matter more transparent for 



protons 



That part of nuclear shadowing which comes from diffractive excitation of the valence 
quark component of the projectile hadron corresponds in terms of the triple-Regge phe- 
nomenology to the PPR term in the diffractive cross section. In eigenstate representation 
for the interaction Hamiltonian the same effect comes from the dependence of the elastic 
amplitude on positions of the valence quarks in the impact parameter plane 0. On top 
of that, the projectile hadron can dissociate via gluon radiation which corresponds to the 
triple-Pomeron term in diffraction. It can also be interpreted in the infinite momentum 
frame of the nucleus as a reduction of the density of gluons which interact with the hadron. 
This relation gives a hint how to approach the problem of gluon shadowing at small x for 
soft gluons. 

Let us model this situation in eigenstate representation with two Fock states for the 
projectile hadron, 

\h) = {l-w)\h), + w\h)G , (114) 

where \h)y and |/i)g are the components without (only valence quarks) and with gluons 
which can be resolved at the soft scale. We assume them to be eigenstates of interaction 
with eigenvalues and ac respectively. The relative weights are controlled by the parameter 
w. The hadron-nucleon and hadron-nucleus total cross sections can be represented as fT^, [5(1, 



(^tot=^v + wAa , (115) 



where Aa = ac — cTv, and 



1 - exp ( a^T{b) 



+ w 



exp ( - ^ (Tv T{b)j - exp ^ ac T{h) 



16) 



This cross section is smaller than one given by the eikonal Glauber approximation 
and the difference is known as Gribov's inelastic corrections IBBI. The Glauber's cross 



section can be corrected by replacing the nuclear thickness function by a reduced one. 
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T(6) =^ T{h) < T{b), which is related to the reduced gluon density in the nucleus, 



GAix,b) T{b) 



;ii7) 



GN{x,b) T{b) ■ 

Thus, nuclear shadowing for soft gluons can be evaluated comparing the total cross 



section (|116|) with the modified Glauber approximation. 
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Expanding both expressions in small parameters AaT and dt, AT, where AT{b) = T{b) 
T{b), (they are indeed small even for heavy nuclei) we get. 
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We left here only the leading terms and omitted for the sake of simplicity the terms con- 
taining higher powers of w. 

According to relation ( |A. 10| ) w(A(j)^/167r is the forward cross section of diffractive 
gluon radiation which corresponds to the triple- Pomeron part of the diffraction cross section 
hN —y XN. Therefore, the correction ( 119 ) to the nuclear thickness function can be 
expressed in terms of the effective cross section. 



eff 



w (Act) 

Otot 



16 n A3p{hN ^ XN) In 



M2 



M2 



(120) 



where M^^^ ^ 2\/?)s/{mM Ra) is the upper cut off for the diffractive mass spectrum im- 
posed by the the nuclear formfactor. The bottom cut off depends on M^-dependence for 
the triple-Pomeron diffraction at small masses which is poorly known. At high energies 
under consideration this uncertainty related to the choice of M^j„ is quite small. We fix 
M^i^ = 2GeV. 

Within an approximate Regge factorization scheme A^p{hN XN) defined in ( P^ ) is 
an universal constant (see, however, (^)). Therefore, the driving term in ( [L19| ) and gluon 
shadowing are independent of our choice for hadron h, a result which could be expected. 

Data on diffractive reaction pp ^ fix the triple-Pomeron coupling {e.g. see in 
PB| , |B^) with much better certainty than for other reactions (including data for diffractive 
DIS). The value of A^p varies from 0.075 GeV~^ at medium high energies to 0.025 GeV~^ at 
Tevatron energy (see (^Tf)). Correspondingly, the effective cross section for A ^ 200 ranges 
as (Te// ~ 3.5 — 5.5 m6. This is an order of magnitude smaller than the value used in [^] at 
high Q^. It is very improbable that cxe// can grow (so much!) with Q^. 

It is silently assumed in Eq. ( p.l6| ) that the energy is sufficiently high to freeze the 
fluctuations, i.e. there is no mixing between the Fock components during propagation 
through the nucleus. If, however, the energy is not high, or the effective mass of the 
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excitation is too large, one should take care of interferences and represent ( p.17] ), (|1 19|) in 
the form (compare to [70(1 ) 



Ga{x] 



AGn{x) 
X pa(&, 2i)pa(&, 22)exp 



1 - Stt Asp{pp pX) Re J d% J -j^ J dzi J dz2 6(^2 - zi 



Mi 



22 



exp 



(121) 



where aabs = ^<^i and we exponentiated the expression in square brackets in the r.h.s. of 
The important difference between (|121|) and the usual expression 



69, ffOl for inelastic 



corrections is absence of absorption for the initial [z < Zi) and final {z > Z2) protons in 
(|121|) . This is a natural result, since proton absorption (mostly of eikonal type) has no 
relevance to gluon shadowing. 



Absorption aabs = ^cr in intermediate state {zi > z > Z2) is much smaller than a 



NN 
tot 

and is related to the amplitude of diffractive gluon radiation (see ( |A.1CI| )). One can estimate 
aabs assuming Regge factorization. In this case Ue// is universal and can be applied even 
to a quark, i.e. h = q. This makes sense in our model due to short range nature of the 
nonperturbative gluon interactions. 

demanding Eq. (|121|) to reproduce correctly the "frozen" limit of 0. This needs 

o'abs = (^eff, as was actually guessed in p][] . 

However, the discussion following Eq. (|4^) shows that after it is averaged over the quark- 
gluon separation the absorptive cross section gains an extra factor, aabs = 1-5 CTg//. 

We performed numerical estimates for A = 200, 64 and 12 assuming a constant nuclear 
density Pa('^) = Po ©(-Ra — 1^) with po = 0.157/m~^ and Ra = 1.15 A^^^ fm. In this case 
integration in ( |121| ) can be performed analytically and the result reads. 
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(122) 



where v = ^aejf po Ra, 7 = 0.5772 is the Euler constant, and Ei{z) is the integral expo- 
nential function. The value of x can be evaluated as x = 4 {k'^)/s, where {k"^) ~ l/fog is the 
mean transverse momentum squared in the quark - gluon system. 



The results of numerical calculations with Eq. |122| for gluon shadowing are depicted in 
Fig. § by thin solid curves for lead, copper and carbon (from bottom to top) as function 
of X. Shadowing for soft gluons turns out to be much weaker than predicted in ^ for 
high Q^. This contradicts the natural expectation that the softer gluons are, the stronger 
shadowing should be. 
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X 

Figure 6: The same as in Fig. ^ but for soft gluons. The thin curves 
are obtained with ( \12^) using data for the triple-Pomeron contribution 
to diffraction pp pX . The thick curves are predicted using the Green 
function method. 



(a) The Green function formahsm 

One can also use the Green function formalism to calculate nuclear shadowing for soft 
gluon radiation. It provides a better treatment of multiple interactions and phase shifts in 
intermediate state. In contrast to the above approach which uses a constant average value 
for (Jeff, in the Green function formalism the absorption cross section as well as the phase 
shift are functions of longitudinal coordinate. This is also a parameter-free description, all 
the unknowns have already been fixed by comparison with data. 

As usual, we treat shadowing for soft gluons as a contribution of the gluonic Fock com- 
ponent to shadowing of the projectile-nucleus total cross section. One can use as a soft 
projectile a real photon, a meson, even a single quark. Indeed, the mean quark-gluon sepa- 
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ration 1/bo ~ 0.3 /m is much smaller that the quark-antiquark separation in a light meson 
or a qq fluctuation of a photon. For this reason one can neglect in (^9|) the interference 
between the amplitudes of gluon radiation by the q and q. Since the gluon contribution to 
the cross section corresponds to the difference between total cross sections for \qqG) and \qq) 
components, the quark spectator cancels out and the radiation cross section is controlled 
by the quark-gluon wave function and color octet (GG) dipole cross section. 

Thus, the contribution to the total hadron-nucleus cross section which comes from gluon 
radiation has the form, 

hA } daa 



^G= ^ d'bPiaG,b) , (123) 

X 

where 

oo 

P{aG,b) = J dzpAib,z) J (fr '^qGir^ac)'^ (TGGij.s) 

— oo 

1 °° 

- 2^^/ dzidz2Q{z2- zi) pAib,Zi) pAib,Z2) J cfri(fr2 (124) 

— oo 

X ^*G(f2, aG)crGGir2, s) GGG{r2, Z2\ ri, Zx) OGG^rx, s) ^gG(n, oig) . 

Here the energy and Bjorken x are related as s = ^m^qv = Ab^/x. The explicit solution for 
the Green function GGG{r2, Z2; ri, zi) in the case of (JGG{r,s) = CGG{s)r'^ and a constant 
nuclear density is given by Eq. (p5|). Note that the approximation for the dipole cross 
section is justified by the small value of (r^) = 1/6q ^ 0.1 /m^. 
Integrations in ( |124| ) can be performed analytically, 

P{aG,b) = ^Reln{W) , (125) 



where 



A2 I h-2 

W = ch{n Az) + ° sh{n Az) , (126) 



Az = 1 JR\ - 62 . (127) 



We use here the same notations as in Eqs. (|109|) - (|110|) . 



The results of calculations are depicted in Fig. | by thick curves for lead, copper and 
carbon (from bottom to top). They demonstrate about the same magnitude of shadowing 
as was calculated above using hadronic basis. However, the onset of shadowing is delayed 
down to X < 0.01. We believe that this result is trustable since the Green function approach 
treats phase shifts and attenuation in nuclear matter more consistently. 

Comparing predicted shadowing for soft gluons in Fig. ^ and one at = AGeV in 
Fig. ^ we arrive at a surprising conclusion that shadowing is independent of scale. A small 
difference is within the accuracy of calculations. This is a nontrivial result since calculations 
were done using very different approximations. Shadowing of hard gluons was estimated 



34 



assuming that the qq pair is squeezed to a size ~ l/Q much smaller than the transverse 
separation between the gluon and the qq. On the contrary, radiation of soft gluons is 
dominated by configurations with a distant q and q surrounded by small gluon clouds. The 
fact that shadowing appears to be the same is a result of existence of the semihard scale 
(which should be compared with Qlfj < <5^/4). At larger virtualities shadowing decreases 
as one can see from comparison of = 4 GeV^ with 16 GeV^ in Fig. ^ 



3.3 Nonperturbative effects in the transverse momentum distri- 
bution of gluon bremsstrahlung 

As soon as the strength of the nonperturbative quark-gluon interaction is fixed, we are in a 
position to calculate the cross section of gluon bremsstrahlung for a high energy quark inter- 
acting with a nucleon or a nuclear target and to compare the results with the perturbative 
QCD calculations p6[] . 



3.3.1 Nucleon target 



The transverse momentum distribution of soft gluons (a^ <^ 1) reads p6| , 
da 1 



2(27r) 



d rid r2 ^^^(^1' "g) ^gG(^2, "g) exp i krin - ^2 



(128) 



Here the overline means that we sum over all possible polarizations of the radiated gluon 
and recoil quark and average over the polarization of the initial quark. In our model for the 
quark-gluon distribution function including nonperturbative effects we get. 
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(129) 



The cross section o"gg('") in ( |128|) has the form 
We performed calculations for the transverse momentum distribution of gluons for two 
parameterizations of the dipole cross section, 

(I) one which is given by (|30D which is constant at ^ Pq. For the sake of convenience we 
change the notation here, = 2/pq = 0.125 GeV~'^; 

(II) the dipole approximation (|) with C = ao s^/2. Only this parameterization is used for 
nuclear targets because it allows to perform integrations analytically (of course one can do 
numerical calculation for any shape of the cross section). 

Correspondingly, we obtain for the differential radiation cross section. 
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and Ei(z) is the exponential integral function. 

In the case of parameterization II it is convenient to represent the dipole cross section 
in the form, 
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Then the differential cross section reads, 
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The results of calculations for variants I and II are depicted in Fig. |^ by solid and dashed 
curves respectively. The two upper curves correspond to perturbative calculations, while the 
two bottom ones include the nonperturbative effects. The strong interaction between gluon 
and quark leads to a substantial decrease in the mean transverse size of the quark-gluon 
fluctuation. Therefore, the mean transverse momentum of the radiated gluons increases. 
The nonperturbative interaction has especially strong effect at small transverse momentum 
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Figure 7: Transverse momentum distribution for gluon bremsstrahlung 
by a quark scattering on a nucleon target. The solid and dashed curves 
correspond to parameterizations I and II for the dipole cross section, 
respectively. The upper curves show the results of the perturbative QCD 
predictions the bottom curves correspond to the full calculation 
including the nonperturbative interaction of the radiated gluon. 



kx, where the radiation cross section turns out to be suppressed by almost two orders of 
magnitude compared to the perturbative QCD expectations. 

Note that intensive gluon radiation originating from multiple nucleon interactions in rel- 
ativistic heavy ion collisions is found 1 78, 79 1 to be an important alternative source for sup- 
pression of charmonium production rate and is able to explain the corresponding data from 
the NA50 experiment at CERN SPS. The found strong suppression of gluon bremsstrahlung 
by the nonperturbative interaction relevant only to small a <^ 1. However, it may substan- 
tially reduce the influence of prompt gluons on charmonium production if is important at 
large a as well. This is to be checked. 
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3.3.2 Nuclear targets 



In the case of nuclear targets Eq. (|128|) holds, but aGcir) has the form, 



a^^(r) = 2 J d^B S^l- exp 



-la,,ir)TiB) 



(135) 



Our calculations for gluon radiation in the interaction of a quark with a nuclear target 
are performed only in the parameterization II for the sake of simplicity. For heavy nuclei 
this approximation can be quite good due to a strong color filtering effect which diminishes 
the contribution from large size dipoles. The transverse momentum distribution has the 
form. 



d afj 8 as 



where 



J d^BF{kl,blS\B)) , (136) 



S\B) = laos'TiB) . (137) 
o 

For numerical calculations we use the approximation of constant nuclear density, pA{f) = 
3A/ {4:7rR\) Q{Ra — r). The results for the radiation cross section per bound nucleon with 
(solid curve) and without (dashed) the nonperturbative effects are compared in Figs. || and 
^ for copper and lead targets respectively. Obviously the nonperturbative interaction 
generates very large nuclear effects. 

The nuclear effects are emphasized by a direct comparison in Figs. 10 and for different 



targets, a nucleon, copper and lead, including and excluding the nonperturbative interaction 
respectively. We see that the difference between a free and a bound nucleon at small 
fcy is substantially reduced by the nonperturbative interaction. Indeed, the interaction 
squeezes the quark-gluon fluctuation and reduces the nuclear effects. Besides, the region of 
antishadowing is pushed to larger values of kr- 

This manifestation of the nonperturbative interaction implies that gluon saturation 
which is an ultimate form of shadowing should happen with a smaller gluon density com- 



pared to the expectations El, p3] based on perturbative calculations. On the other hand. 



the saturation region spreads up to higher values of kx- 



4 Summary and outlook 

We explicitly introduced a nonperturbative interaction between partons into the evolution 
equation for the Green function of a system of quarks and gluons. The shape of the qq 
potential is chosen to reproduce the light cone wave function of mesons. The magnitude of 
the potential is adjusted to reproduce data for photoapsorptive cross sections on nucleons 
and nuclei and data on diffractive dissociation of photons into qq pairs. 

Based on theoretical arguments and experimental facts we expect a much stronger in- 
teraction for a quark-gluon pair than for a quark-antiquark pair. Indeed, data on diffractive 
dissociation of hadrons and photons into high mass states show that the cross section is 
amazingly small, what is usually phrased as evidence that the triple-Pomeron coupling is 
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Figure 8: The differential cross section per bound nucleon of soft gluon 
brenisstralilung in quark-copper collisions. The solid and dashed curves 
correspond to calculations with and without the nonperturbative effects 
respectively. 



small. We have performed calculations for diffractive gluon radiation (responsible for the 
production of high mass excitations) including the nonperturbative effects, and fixed the 
strength of the quark-gluon potential. We found a very simple and intuitive way to get the 
same results as direct calculations of Feynman diagrams. Both approaches lead to the same 
diffractive cross section which in the limit of perturbative QCD coincides with the result 
of a recent calculation |19] for the process 7*A^ qqGN. We conclude that the previous 
analogous calculations |]18| are incorrect. 

We adjusted the quark-gluon potential to data for the diffractive reaction pp Xp which 
have the best accuracy and cover the largest range of energies and masses. We predicted the 
single diffractive cross sections for pions, kaons and photons and find a substantial violation 
of Regge factorization. 

We calculated nuclear shadowing for longitudinally polarized photons which are known 
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Figure 9: The same as in Fig. 13, but for a lead target. 



to serve as a sensitive probe for the gluon distribution, using the Green function technique 
developed in 



12 



describing the evolution of a qqG system propagating through nuclear 
matter. The evolution equation includes the phase shift which depends on the effective 
mass of the fluctuation, nuclear attenuation which depends on the transverse separation 
and energy, and the distribution over transverse separation and longitudinal momenta of 
the partons which is essentially affected by the nonperturbative interaction of the gluon. 
The latter substantially reduces the effect of nuclear shadowing of gluons. We have found 
an X dependence for gluons which is quite different from that for quarks. These differences 
are far beyond the simple Casimir factor 9/4. 

Nuclear shadowing for soft gluons is essentially controlled by the nonperturbative effects. 
It turns out to be rather weak similar to what is found at ~ AGeV"^. Such a scale in- 
variance at low and medium high virtualities is a consequence of the strong nonperturbative 
interaction of gluons which introduces a semihard scale ~ 4 6q = 1.7 GeV"^. 

The nonperturbative interaction changes dramatically the transverse momentum distri- 
bution of gluon bremsstrahlung by a high energy quark interacting with a nucleon or a 
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Figure 10: Comparison of the cross sections of gluon radiation per nu- 
cleon in the perturbative QCD hmit for coUisions of a quark with a nu- 
cleon (sohd curve), copper (dashed curve) and lead (dotted curve) versus 
the transverse momentum squared of the gluon. 



nucleus. The gluon radiation cross section at small kx turns out to be suppressed by nearly 
two orders of magnitude compared to the expectations from perturbative QCD ||5^, 
Although these results concern the gluons radiated with ac — > 0, it might also suppress 
gluon bremsstrahlung at larger ac which is predicted ||78| to contribute to the break up of 
charmonia produced in relativistic heavy ion collisions. 

This effect is especially strong for nuclear targets where the nonperturbative interaction 
of radiated gluons creates a forward minimum in the transverse momentum distribution. 
This suppression is an additional contribution to nuclear shadowing calculated perturba- 
tively in |^0|, |36| which also leads to a suppression of small transverse momenta. The results 
of our calculations presented in Figs. |^ include both phenomena. 

Nuclear shadowing for small transverse momenta of the radiated gluons is the same effect 
as the saturation of parton densities at small x in nuclei as seen in the infinite momentum 
frame of the nucleus. This phenomenon is expected to be extremely important for the 
problem of quark-gluon plasma formation in relativistic heavy ion collisions. On the one 
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Figure 11: The same as in Fig. [I^, but the nonperturbative interaction 
of gluons is included. 



hand, a growth of the mean transverse momentum of radiated gluons increases the produced 
transverse energy, on the other hand, it leads to a higher probability for such gluons to escape 



the interaction region without collisions, i.e. the gluon gas may not reach equilibrium |80 
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Appendix A Diffraction 

A.l General consideration 

In this section we present a general analysis of diffraction based on the eigenstate decom- 
position. 

The off-diagonal diffractive scattering is a direct consequence of the fact that the inter- 



acting particles (hadrons, photon) are not eigenstates of the interaction Hamiltonian [^T], ^ 
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They can be decomposed in a complete set of such eigenstates \k) 

\h) = Y.Ct\k). 



(A.l) 



where are the amphtudes for the decomposition which obey the orthogonahty conditions, 



h 



h' 
5lk ■ 



(A.2) 



We denote by fk = i o'k/'^ the eigenvalues of the elastic amplitude operator /. We assume 
that the amplitude is integrated over impact parameter, i.e. that the forward scattering 
elastic amplitude is normalized as = 4:7idak/dt\t=o- We can then express the hadronic 
amplitudes, the elastic feiihh) and off diagonal diffractive fddihh') amplitudes as. 



Uhh) = 2tY: \ci 



CTfc 



2t (a) 



(A.3) 



(A.4) 



Note that if all the eigen amplitudes are equal the diffractive amplitude ( A.4 ) vanishes due to 
the orthogonality relation, ([A .21) . The physical reason is obvious. If all the fk are identical 
the interaction does not affect the coherence between the different eigen components \k) 
of the projectile hadron \h). Therefore, off diagonal transitions are possible only due to 
differences between the f^s. For instance, in the two channel case. 



fad{hh') = 2z{C'^yC^iai-a2) . 



(A.5) 



If one sums over all final states in the diffractive cross section one can use the complete- 
ness condition ( |A.2|) . Excluding the elastic channels one gets 0, |^, 



IGtt 



da. 



dd 



dt 



t=o 



0", 



(A.6) 



This formula is valid only for the total (forward) diffractive cross section and cannot be 
used for exclusive channels. 



A.2 Diffractive excitation of a quark, q — > qG 



In this case we can restrict ourselves to the first two Fock components of the quark, a 
bare quark \q) and \qG). Therefore, we can use Eq. ( A.5 ). Thus, we arrive at the following 
expression for the forward amplitude of diffractive dissociation into a qG pair with transverse 
separation p, 



fdd{q qG) 



gT=0 



a. 



qG 



(P) 



(A.7) 
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Both cross sections, and a, are infra-red divergent, but this divergence is obviously the 
same and cancels in ( [A.7| ). 



To regulate the divergence we can introduce a small gluon mass m^, which will not 
enter the final result, and impose that for separations r ^ l/mc the dipole cross section 
is given by the additive quark limit, <Jqq{r ^ l/mc) = 2aq. To find the convergent part of 
(^qcip) — o'q "we can make use of Eq. (W^). Let us choose in (W^) ri <ti l/mc and r2 ^ l/mc- 
Then the l.h.s. of { ^7\) saturates at cr^ + (Jqciri). Here CTqciri) is different from ag due to 
the color dipole moment of the qG system, i.e. due to ri 7^ 0. Then (^Tj) is modified to, 

9 r . . „ 1 2 



9 f 12 

(Tg + (Tqcin) = - ^crgg{ri) + 2(Tg| - - ag . (A.8) 



From this relation we obtain the combination of cross sections at the r.h.s. of Eq. ( |A.7| ) 
which takes the form, 



fddiq qG) 



= i^gG{a,p)lagg{p) . (A.9) 

qT=0 O 



Thus, we derived Eq. (^Oj) in a simple and intuitive way. A more formal derivation based 
on direct calculation of Feynman diagrams is presented in Appendix B.l. 



A. 3 DifFractive gluon radiation by a qq pair 



The diffractive amplitude of gluon radiation by a qq pair, qq qqG, can be easily derived 
in this approach. We restrict ourselves to two Fock components \qq) and \qqG). Then the 
distribution amplitudes get the meaning of distribution functions for these Fock states, 
namely \E'gg(ri — ^2, a) and "^cqqipi, P2, «, «g), where the transverse coordinates are defined 
in (0). Summation over k in ( [A.l| ) - ( [A.4| ) now means integration over the transverse 
separations and summation over the Fock components. According to ( |A.4|) - ( |A.5| ) the 
diffractive amplitude fddiQI QqG) reads, 

fddiqq qqG) = 2 J d^pi £p2 "^qqcipi, p2, a, ac) [acqqipi, P2) - (^qqiPi - P2)] ■ (A. 10) 

Here we make use of the obvious relation C^'^{f) = S{f). The total cross sections for the two 
Fock components qq and Gqq are introduced in (^) and (P?]). 

The distribution amplitude for the Gqq fluctuation in the limit of ac — > is easily 
guessed. Indeed, in this limit the impact parameters of the q and q are not affected by gluon 
radiation. Therefore, the "^qqc should be a product of the qq distribution function in the 
projectile hadron (photon) times the sum of the gluon distribution amplitudes corresponding 
to radiation of the gluon by q or q, 



'^qqci.Pl, P2, a, ac) = ^qq{,Pl " P2, «) 



a 



(A.ll) 



where '^gg and ^gc are defined in (|T8D, (|^) and in (^) respectively. Thus, we have arrived 
at Eq. (^). A more formal derivation based on the calculation of Feynman graphs is 
presented in the next Appendix. 
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After integration over (pi + P2) in ( A.lOj ) the amplitude of diffractive gluon radiation 
turns out to be proportional to the difference Aa^pi — P2) between the cross sections of the 
colorless systems Gqq and qq. This is a straightforward consequence of the general property 
of off-diagonal diffractive amplitudes given in ( |A.5| ). 

These conclusions are also valid for diffractive gluon radiation by a photon 7 — >■ 
qqG N . At first glance presence of a third channel, the photon, may change the situation 
and gluon radiation amplitude may not be proportional to Ao". This is not true, however, 
since the relative weights of the qq and qqG components of the photon are the same as above 
as soon as they are generated perturbatively. 

In the limit of purely perturbative interactions the same result as our Eq. (l83|) was 
obtained recently in (Eq. (3.4)). However, the cross section for diffractive gluon radiation 
derived earlier in ||18| (Eq. (60)) is not proportional to (Aa)^, but contains a linear term. 
We think that this is a consequence of improper application of relation ( |A.(j| ) to an exclusive 
channel. 



A. 4 Diffractive electromagnetic radiation 

The forward amplitude for photon (real or virtual) radiation by a quark is similar to that 
for gluon radiation ( |A.7|) , except that the photon does not interact strongly and one has to 
replace by cTg, 



fdd{q qi* 



gT=0 







(A.12) 



Thus, in order to radiate a photon the quark has to get a kick from the target, no radiation 
happens if the momentum transfer to the target is zero. 

This conclusion is different from the expectation for diffractive Drell-Yan pair production 
of |2^. The latter was based on the conventional formula ( |A.l(j| ) which cannot be used for 
an exclusive channel (as well as for gluon radiation). Therefore, the diffractive Drell-Yan 
cross section should be much smaller than estimated in |^ . 

Nevertheless, a hadron as a whole can radiate diffractively a photon without momentum 
transfer as two of its quarks can participate in diffractive scattering, each of them may 
getting a momentum transfer, while the total momentum transfer is zero. 



Appendix B Diffraction: Feynman diagrams 
B.l qN^qGN 

For the example of diffractive excitation of a quark, 

qN-^qGN, (B.l) 

we demonstrate in the following the techniques and approximations we use for the calculation 
of more complicated diffractive processes. 
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We use the following notations for the kinematics of kx and pt are the transverse 

momenta of the final gluon and quark respectively; a is the fraction of the initial light- 
cone momentum carried by the gluon; = kx + Pt is the total transverse momentum 
of the final quark and gluon, and kt = (1 — ajkx — otpr appears further on, when the 
transverse separations re = 6 + (1 — a)p and = 6 + ap are inserted: kx ■ tg + Pt ■ rq = 

{kr + Pt) ■ 6 + ((1 - a)kT + ^Pt) ■ P- 

We normalize the amplitude of ( |B.1[ ) according to 

dm a cI^kt d^qT -J ^~^^s 



i ^ Tr \aI {qr, kt, «) (gr, ^t, a)] , (B.2) 



where 



/l^"^ = (g^)^4g% (B.3) 

and g^*^^^ are the color spinors of the quark in the initial and final states; s is the color index 
of the radiated gluon. 

We assume that at high energies one can neglect the ratio of the real to imaginary parts 
of the amplitude for reaction (p.l| ). Then one can apply the generalized optical theorem 
{Cutkosky rules |p8|), 

A{a^b) = -^A^{b^c)A{a^c) . (B.4) 
2 c 

here includes not only a sum over intermediate channels, but also an integration over the 

c 

intermediate particle momenta. 

To simplify the problem we switch to the impact parameter representation, 

A{b,p) = ^^^^^ jd'^qd'^KTA{qT,KT)exp(^-iqTb-iKTp) ■ (B.5) 

Since the initial impact parameters are preserved during the interaction we sum only over 
intermediate channels in this representation. 

We use the Born approximation, i.e. the lowest order in a^, for the sake of clarity, and 
generalization is straightforward. In this case and for a = {qN}, b = {qGN} only two 
intermediate states are possible in ( p.4| ): ci = {qN^} and C2 = {qG N^}, where iVg = |3g)8 
is the octet color state of the 3q system produced when the nucleon absorbs the exchanged 
gluon. 

One should sum in (|B.4|) over all excitations / of the N^, 

MqN ^ qGN) = ^ [A^qGN ^ qN*,) i(giV ^ giV*) 



/ 



+ E A^qGN ^ qGN;) A^qN ^ qGN*,) 



(B.6) 
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Here s' is the color index of the gluon in the intermediate state. 

We skip the simple but lengthy details of calculation of the amplitudes on the r.h.s of 



(|B.6|) and present only the results. 

A(qN ^ qN;) = Tr (f 



(B.7) 



A',{qN^qGN; 



(B.8) 



A,{qGN qN; 



rsp '^p 



f Irih 



AsAqGN qCN;) = S,,,Tr If l)+lfss'r (f lr{h) 



(B.9) 

Here hi = b; b2 = b — a p are the impact parameters of the projectile and ejectile quarks 
in reaction (p.l| ), respectively; 63 = + (1 — a) p is the impact parameter of the radiated 
gluon; p is the transverse separation inside the qG system, and b is the distance from its 
center of gravity to the nucleon target; \l/gG(«, p) is the distribution function for the qG pair; 
Aj. = 2 Tr are the Gell-Mann matrices; frsp is the structure constant for the SU{3) group. 
The matrices 7r(&fc), {k = 1,2,3) are the operators in coordinate and color space for the 
target quarks. 



%{bk 



X0) = - d'q 



(B.ll) 



(B.12) 



g2 + A2 

where is the transverse distance between the j-th valence quark of the target nucleon 
and its center of gravity; the matrices rj^^^ act on the color indices of this quark. The 
matrix elements (/|7r(&A;)K) between the initial i = N and final / = states are expressed 
through the wave functions of these states. The effective infra-red cut off A in (|B.12|) does 
not affect our results, which are infra-red stable due to color screening effects. 
Substitution of (|R7D-( plO|) into ( TO results in. 



MqN ^ qGN) 



= 2 V^'^^^r' \^rr'ibi,bi) - V |$rr'(&2,&l) + i frspTpTri |$„/(63,6i) 

+ TrTsTr> \(^Tr'Q^lM) - V |$rr'(^2,&2) -ifrspTpTr' |$rr'(^2,&3) 

- i fss'r Ts> V |$rr'(^l, ^s) " V Ts' |$rr'(^2, h) - i fs'r'p Tp ^rr'ih, h) } 

73 



X 



*<?G(a,p) 



(B.13) 
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where 



f 



(B.14) 



We sum in ( |B.13|) over all excitations of the two color octet states of the 3q system. To 
have a complete set of states we have to include also color singlet and decuplet |3g) states. 
As these states cannot be produced via single gluon exchange, they do not contribute and 
we can simply extend the summation in (|B.14| ) to the complete set of states and get, 



lr{hk)lr'{hl 



(B.15) 



In the matrix element ( |B.15| ) we average over color indices of the valence quarks and 
their relative coordinates in the target nucleon. To do so one should use the relation. 



|3g>i 



Then, Eq. ( [B.15| ) can be represented as. 



6rr' S{bk, hi 



(B.16) 



(B.17) 



where 5(6^, h) is a scalar function of two vector variables. 



S{hM) = ^ I d{s} 



3 ^ ^ 1 ^ 

Xih - Sj) xik - Sj) - - II xih - Sj) xik - Sj. 



2 ^, 



This function is directly related to the qq dipole cross section 



I<f3.({4)r • 

(B.18) 



(T5g(pi -P2)= j d'b [S{b + p,,b + p,) + S{b + P2,b + p2) - 2 S(6 + Pi, 6 + p2)J . (B.19) 

According to ( [B.17|) and ( B.18|) the function $(6^, bi) is symmetric under the replacement 
bk ^ bi. Therefore, the terms proportional to $(61,62) and $(62,61) in ( |B.13D cancel, as 



well as the terms proportional to $(61,63) and $(63,61) 
proportional to $(62, 63) and $(63, 62) add up. 
Making use of the relations. 



At the same time, the terms 



fss'r fs'rp 
frsp Tr 



4/3 
36 



sp 



(B.20) 



we arrive at the final result for the amplitude of diffractive dissociation of a quark {qN 
q G N) in impact parameter representation, 

J /T - \ i3\^ ^ , r 4 r 
As{b,p,a) = TsWgG{a,p) 



16 



5(61,60-5(62,62) 



+ 3 



5(61, 63) -5(63, 63)]} 



(B.21) 



48 



The diffraction amplitude in momentum representation reads, 

^sigr^f^TjOi) = j^—j^ J d%d^pAs{b,p,a)exp (iqTb + iKTpj • (B.22) 



Using ( |B.19D and the above mentioned symmetry of S{pi,p2) we obtain a very simple 
expression for the forward {qx = 0) diffraction amplitude which is related to the dipole cross 
section, 

^9^3 



As{0,KT,a) 



Ts / d p-^qG{.a,p)aqq{p)e 



32 (27r)2 

Eventually, the forward diffractive dissociation cross section of a quark reads, 
da 

~ 3 



dilna) d^qx 



ijt=0 



^ j d^KT E TriI(0,/?T,a)i.(0,/?T, 



a] 



1 



(4vr) 



d^p 



9 



We should emphasize that all above calculations are done for an arbitrary a. 



(B.23) 



(B.24) 



B.2 Diffractive gluon radiation by a qq pair 

Gluon radiation is an important contribution to the diffractive dissociation of a (virtual) 
photon, 

-f* N qqGN . (B.l) 
In analogy to the previous section we make use of the generalized unitarity relation, 

MYN ^ qqGN) = ^ E [MmGN ^ ggiV*) i(7*iV ^ ggiV*) 

^ / 

+ J2^sAqqGN^ qqGN;) AA^*N-^ qqGN;) , (B.2) 

s' 

where the amplitudes are defined as follows, 

A{j*N qqN;) = k (/ Uh) ^) + (/ Uh) ^)] ^5,(pi - P2, a) \qq), ; (B.3) 



AAl'N^qqN; 



iV3 



fs'rp {rp [{f\lr{bi)\i) - {f\lr{h)\i)] *</g(pi) 



f\lr{b2)\i) - {f\lr{h)\i)] ^yG(P2)} ^g-g(Pl -P2,tt) \qq) i 



(B.4) 



AiqqGN^qqN;) 



fsrp {rp [{f\lr{bi)\i) - {f\lr{h)\i)] ^gG(Pl) 

/|7r(&2)|^)-(/|7r(&3)|^)] ^ qG{p2)} ] 
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(B.5) 



Ass'{qqGN -> qqGN*) = (^f %{bi) 5ss' + (/ 7r(&2) 



(B.6) 



Here bi = b + fi, b2 = b + f2, b-^ = b + p are the impact parameters of the quark, antiquark 
and gluon respectively; b is the photon impact parameter; pi,2 = f> — n,2; ^qq and 
are spatial and color parts of the gg-component of the photon wave function, respectively. 
The matrices = Ar/2 and tv = A*/2 act on the color indices of quark and antiquark 
respectively. The indices s, s' mark the color states of the gluons in intermediate and final 
states. 

Note that the condition of color neutrality of the singlet state \qq) leads to the relation. 



{Tr + Tr) \qq) = 



(B.7) 



Substitution of (|B.2|) - (|B.6|) into (|B.1|) leads to the following expression for the amplitude 
of diffractive dissociation of the photon. 



A{-f*N ^ qqGN) 



16 



{ fsrp [rp Tr + TrTp] S (61 , 61 ) - S (63 , 61 ) 



+ i fsrp [Tp Tr + Tr Tp] s{bi, 62) - s(&3, ^2) + fss'r fs'rp Tp 3(63, 61) - 5(63, 63) ) ^gG(Pl) 
+ (i fsrp [Tpfr + frfp] 3(62, ^2) " -5(62, ^s) + i fsr ■p [Tp Tr + Tr Tp 

] S{b2,b,) - S{b,,b,) 

+ fss'r fs'rpfp 3(63,62) - 5(63,63) ) ^gG(P2)} *g-g(Pl - P2, a) , (B.8) 

where we made use of the completeness condition, J2f I/) (/I = 1 (see Appendix B.l). 

In order to simplify Eq. (p.8| ) we apply a few relations as follows. Since fsrp = —fspr we 
find 

fsrp \Tp Tr ~l~ Tr Tp] fsrp \_Tp Tr ~l~ Tr Tp] . (B.9) 

Then, relying on the condition ( p.7|) we find. 



[TpTr + TrTp) \qq) ^ 



2TpTr \qq)-^ 
-2TpTr \qq)^ ; 



(B.IO) 



{TpTr + TrTp) \qq)^ = 2TpTr\qq)^ 

= -2fpfr \qq)^ . (B.ll) 



We also use the relations 



^ f srp Tp Tr ^ i 

•f - - - ■ 

^ I srp Tp Tr 2 ' 

fss'r fs'rp 3 (J^p , (B.12) 
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and the symmetry condition, s{bk,bi) = s{bi,bk), and eventually arrive at a modified form 
of Eq. 



M-f*N ^ qqGN) 



9^3 
16 



[ts ^gG(Pi) + fs ^gG{p2)] \qq)i - P2, 



a] 



X 



(B.13) 



s(&2, bs) + s{bi, 63) - s{bi, 62) - 5(63, &3 
The last factor in square brackets can be represented as, 

P{h,b2]h) = S{b2,b3) + S{bi,h) - S{bub2) - S{h,h) 

= ^ {[s{h,h) + S{b2,h) ~ 2 S{h,h)] - [5(62, &2) + 5(63, &3)- 2 5(62, &3) 



(B.14) 



Then, the forward diffraction amplitude {qt = 0) in impact parameter representation has 
the form, 



1 

2^ 



(fb As{b, P1P2) 



4 TT 



— r,^,G(pi) + ^r,^,G(p2) 



199)1 ^g<?(Pi -P2,«) , (B.15) 



where T,{pi,p2) is introduced in (|73|). 

From ( [B.lSp one easily gets the forward diffractive cross section. 



da{YN qqGN) 



dilnac) dqr 



qr = 



d^Pi d'^p2 da 



X 



(4 7r)2 
*,g(pi) - ^9g(P2) 



*w(Pl - P2,tt) 



S(P1,P2; 



(B.16) 



B.3 Diffractive photon radiation, qN — ^ 79 AT 



Diffractive electromagnetic radiation is calculated in analogy to what was done in Appendix 
B.l for gluon radiation. Since the photon does not interact with the gluonic field of the 
target the structure of all the amphtudes in the relation, 



A{qN-^qjN) = 1 ^ [it(^^;v ^ giV*) i(giV giV* 
+ A\q^N q^N;) i(giV q^N*,] 

turns out to me much simpler. 

A(qN ^ qN;) = (/ %ih) 



(B.l) 



(B.2) 



51 



(B.3) 



A{qN ^ q^N;) 



A{q-fN ^ q-fN; 



Tr 



f 



f 



(B.4) 



Here bi = b, b2 = b — a p are the impact parameters of the quark before and after radiation 
of the photon; p is the transverse separation between the quark and photon in the final 
state; and a is the fraction of the quark hght cone momentum carried away by the photon. 

a) is the distribution function for the qj fluctuation of the quark. The initial, 
and final, |/), states of the target, as well as the operators ^{bk) {k = 1, 2) are the same as 
in Appendix B.l. 

After substitution of (^) - (IK^ into (|B3D we get. 



A{qN q-fN) = ^ {r, V [$rr'(&l, ^l) - ^rr'ih, h) + ^rr'ih, h) - $„'(62, ^2)] } • (B.5) 

Here the functions ^rr'{bk,bi) are defined in Appendix B.l Then, the amplitude in impact 
parameter representation reads. 



Mb,. 



s{bi,bi) - s{b2M) 



1 r 



s(6, b) — s{b — a p,b — a p) 



(B.6) 



After Fourier transform to the momentum representation we get for the forward diffractive 
amplitude of photon radiation. 



MQT,KT)\g^=o = j^;^ J d^d^pexp (iqrb + iKTp) A{qT,KT] 



. 



(B.7) 



qT=0 



Thus, the direct calculation of Feynman diagrams confirms our previous conclusion (Ap- 
pendix A. 4) that a quark does not diffractively emit electromagnetic radiation if the mo- 
mentum transfer with the target is zero (as different from the statement in [^). A hadron, 
however, can radiate in forward scattering. 



Appendix C The triple-Pomeron coupling 

In the limit of vanishing quark and gluon masses the quark-gluon wave function (^) retains 
only the second term oc $1 which has the form (|27|). Bi-linear combinations of this wave 
function averaged over final polarizations can be represented as follows, 

1 °° 



00 00 



Mp^,<y)■^M) = ^ J J dte-'"^--''^ . (C.2) 



(27r)^ 
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This together with (|30D and (|79D allows to integrate analytically over the coordinates of 
the quarks and the gluon in (178]). Finally integrating over t and u we arrive at, 



G3p{NN XN) 



a,. 



(4vr) 



9 



where x = 6^(0)pq, z = zn = 2{r^)p/pl, and 



Here 



FAx,z) 



In 



+ 2si In 



2 



a;(a; + 2) 

fSiX + S2 



+ 2 S3S4 



2 In 
2 In 



V XSi 

f SiX + S3S4 



Fi{x,z) - F2{x,z) 



{x + l){x + si) 

x{l + X + Si) 
X 

X + 2S4\ 



1 



Si 

S4 

zFo(x, z) 



1 + z' 
2 

2 + 32' 



S2 



XSi 



1 

1 + 2^' 
4 



-In 



X 



S3 



2 + z' 



S5 



4 + 3z ' 



(C.3) 



(C.4) 



(C.5) 
(C.6) 



where 












i = 1 


^ = 


= 2/3, 


(3 = l/z + 2, 


7 = 5 = x/2 + /3 ; 




i = 2 


9 


= 2, 


(3 = l/z + l, 


7 = 5 = x/2 + /3 ; 




i = 3 


9 = 


-10/3, 


(3 = l/z + l, 


6 = x/2 + 13, 


7 = 5 + 1 


i = 4 


9 


= 1, 


(3 = 1/z, 


7 = 5 = x/2 + /3 ; 




z = 5 


9 -- 


= -4, 


(3 = 1/z, 


5 = x/2 + /?, 


7 = (5 + 1 


z = 6 


9 = 


= 5/3, 


/3 = 1/z + 2, 


7 = 5 = x/2 + /? + 1 ; 




i = 7 


9 


= 2, 


(3 = 1/z, 


5 = x/2 + /?, 


7 = 5 + 2 


i = 8 


9 = 


-- S4/3, 


(3 = 1/z + 1/2, 


7 = 5 = x/2 + /3 ; 




i = 9 


9 = 


-2S5/3, 


(3 = l/z + S5/4, 


5 = x/2 + 13, 


7 = 5 + 1 


^ = 10 


9 = 


2S5/3, 


/3 = 1/2 + 1 + S5/4, 


7 = 5 = x/2 + /3 ; 




i = 11 


9 = 


-S4/3, 


(i = l/z- S4/2 


7 = 5 = a;/2 + /5 + S4 ; 




i = 12 


9 = 


-2S5/3, 


/3 = l/z + l-S5/4, 


7 = 5 = x/2 + /5 + S5/2 ; 




i = 13 


9 = 


-2S5/3, 


/? = 1/z - S5/4, 


5 = x/2 + (3 + S5/2, 


7 = 5 + 1 


i = 14 


9 = 


2S4/3, 


/3 = 1/z, 


5 = x/2 + /3 + 1/2, 


7 = x/2 + /? + 



(C.7) 

The effective triple- Pomeron coupling G^pIMN — * XN) for diffractive dissociation of a 
meson M can be calculated in a similar way assuming a Gaussian shape of the quark wave 
function of the meson, 



1^ 



M^qq 



(Ci 
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where i?^ = 8 {r\^)ch/'i- The triple- Pomeron couphng is smaller by a factor 2/3 (different 
number of valence quarks) and has a form similar to ( |C.3| ), 



G^p{MN XN) 



2 a., 



3(47r)2 V8 



9 



(To 



F^{x,zm)-F^\x,zm 



Ml 



(C.9) 



but zm = R'^/Po 7^ and the functions F^2 different too. The expression for zm) 
results from Fi{x, z) via the replacement S3 ^ 1, S4 ^ S2 and z zm- 

The expression for F2{x, Zm) follows from F2{x, z) after moderate modifications in (|C.7|) : 
^(1) = 1, g{3) = -4, g{6) = 2, all g{{} = for i > 8 and 2 ^ zm- 

In the case of diffractive dissociation of a photon, the calculations are more complicated 
since the spatial distribution of quarks in the photon is very different from a Gaussian. 
Nevertheless, it can be represented as a superposition of Gaussians, 



9 l/a^{a) 

2 _ aemNc E 



27r 



i?2 \nR^ 



Then, the effective coupling Gsp{'~fN — > XA^) takes a form similar to (|G.3|) and 



Gspi^N XN) 
Po] 



6(27r) 



where z^ = la^ia) Pn] ^ and 



FUx,z,) = J^Fliix^v). 



Fi{x,z^ 



F^{x, z^) 



(C.IO) 
I), 

(C.ll) 



(C.12) 
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